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Abstract 

The study of solar neutrinos has given since ever a fundamental contribution both to astroparticle 
and to elementary particle physics, offering an ideal test of solar models and offering at the same time 
rS relevant indications on the fundamental interactions among particles. After reviewing the striking 

results of the last two decades, which were determinant to solve the long standing solar neutrino 
puzzle and refine the Standard Solar Model, we focus our attention on the more recent results in this 
field and on the experiments presently running or planned for the near future. The main focus at 
the moment is to improve the knowledge of the mass and mixing pattern and especially to study in 
detail the lowest energy part of the spectrum, which represents most of solar neutrino spectrum but 
is still a partially unexplored realm. We discuss this research project and the way in which present 
and future experiments could contribute to make the theoretical frawemork more complete and stable, 
understanding the origin of some "anomalies" that seem to emerge from the data and contributing to 
answer some present questions, like the exact mechanism of the vacuum to matter transition and the 
solution of the so called solar metallicity problem. 



14.1 Motivations for the solar neutrino study 

The analysis of neutrinos emitted in the fusion processes inside the Sun is one of most significant 
examples of the relevant role played by the study of neutrino properties in elementary particle physics 
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and astrophysics and in creating a link between these two sectors. The pioneering work in the sixties 
PP had the main goal of understanding better the way in which our star shines and to test solar models. 
But, the surprising result of an apparent deficit in the electron neutrino flux reaching the detector 
marked the raise of the so called solar neutrino puzzle, and opened a whole new field of research that 
has been central in elementary particle physics for many decades. 

The experimental results obtained using different techniques in more than thirty years and the 
parallel theoretical advancements confirmed at the end the validity of Pontecorvo's revolutionary idea 
of neutrino oscillation [2], proving in a crystal clear way that neutrinos are massive and oscillating 
particles. This is one of the first pieces of clear evidence of the need to go beyond the Standard Model 
of electroweak interactions and the attempt to accommodate the experimental results about neutrino 
masses and mixing is a test every theory "beyond the Standard Model" has to pass. Therefore, it is 
clear why these results have had a great impact on elementary particle physics and also on cosmological 
models. At the same time, the possibility of measuring directly at least some components of the solar 
neutrino spectrum and of recovering in an indirect way the value of total solar neutrino flux have been 
fundamental for the progressive refinement of the Standard Solar Model (SSM), which has evolved 
during these years and is now in a general good agreement with the solar neutrino experiments. 

Despite the fundamental steps forward made in the last decades, many questions are still open 
about the real nature and the main properties of neutrinos and the exact mixing mechanism, e.g. are 
neutrinos Major ana or Dirac fermions, determination of mass hierarchy and exact mass values, accurate 
determination of the mixing angles, presence of CP violation. The solar neutrino experiments presently 
running or planned for the future can contribute to solve at least some of these puzzles. The new 
frontier in this field is the study of the low energy part of the solar neutrino spectrum which represents 
the great majority of the spectrum, and is still an almost unexplored realm. Some of the challenges 
ahead are: reducing significantly the indetermination on pep and CNO neutrinos and attaching the 
pp solar neutrino measurement. This would be essential to test the stability and consistency of the 
standard explanation of the oscillation mechanism, confirming or definitely disproving the presence 
of discrepancies between theory and experiments, which has lately stimulated a flourishing of models 



introducing the so called "Non Standard Interactions" (Section 14.6.1 ). Once more, these results would 
be of great interest to improve the knowledge both of elementary particle properties and interactions 
and of the astrophysical models of the Sun. They could help also to discriminate between different 
versions of the solar models, for instance of the so-called "solar abundance problem" , and to deepen the 
comparison with the results coming from other studies of solar properties, e.g. from helioseismology. 
This research project would of course imply a further improvement of the already known detection 



techniques and the introduction of new ones (see, for instance, the section 14.7). Also from this point 
of view, solar neutrino physics will continue to give a stimulating contribution both to elementary 
particle physics and to astrophysics. 

14.2 Brief history and solution of the solar neutrino problem 
14.2.1 From Homestake to Super-Kamiokande 



The first experiment built to detect solar neutrinos took place in the Homestake gold mine in South 
Dakota [Ij. The detector consisted of a large tank containing 615 tons of liquid perchloroethylene, 
chosen because it is rich in chlorine and the experiment operated continuously from 1970 until 1994. 
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Neutrinos were detected via the reaction: 

v e + 37 C1 -> 37 Ar + e~ . (14.1) 

The energy threshold of this reaction, Eth = 814 keV, allowed the detection of 7 Be and 8 B (and a 
small signal from the CNO and pep) but not that of pp neutrinos, because of their low maximal energy 
of 0.42 MeV. The radioactive 37 Ar isotopes decay by electron capture with a Ti/ 2 of about 35 days 
into 37 CI*: 

37 Ar + e~ -»• 37 CT + u e . (14.2) 

Once a month, after bubbling helium through the tank, the 37 Ar atoms were extracted and counted. 
The number of atoms created was only about 5 atoms of 3 Ar per month in 615 tons C2CI4. The 
number of detected neutrinos was about 1/3 lower than expected by the Solar Standard Model. This 
discrepancy is the essence of the Solar Neutrino Problem, which has been for many years an important 
puzzle among physicists. 

There were three possible explanations to the Solar Neutrino Problem. The first one was to consider 
that Homestake could be wrong, i.e. the Homestake detector could be inefficient and, in this case, its 
reactions would not have been cpredicted correctly. After all, to detect a handful of atoms per week 
in more than 600 tons of material is not an easy task. The second one was to consider that the SSM 
was not correct, but as helioseismologjfj] started to provide independent tests of solar models the SSM 
passed all tests. Indeed, non-standard solar models constructed ad-hoc to resolve the Solar Neutrino 
Problem seemed very unlikely when scrutinized under the light of helioseismology. The third one, 
and the strangest hypothesis, was to consider that something happens to the neutrinos while traveling 
from the core of the Sun to the Earth. 

The first real time solar neutrino detector, Kamiokande, was built in Japan in 1982-83 [3]. It 
consisted of a large water Cerenkov detector with a total mass of 3000 tons of pure water. In real 
time neutrino experiments scientists study the bluish light produced by the electrons scattered by an 
impinging neutrino according to the following equation: 

v x + e~ — > v x + e~ . (14.3) 

In the Kamiokande detector light is recorded by 1000 photomultiplier tubes (PMT) and the energy 
threshold of the reaction is Eth = 7-5 MeV; therefore only 8 B and hep neutrinos are detected. At the 
beginning of the '90s a much larger version of the detector was built, Super-Kamiokande, where the 
active mass was 50000 tons of pure water viewed by 11200 PMTs. In Super-Kamiokande the energy 
threshold was lowered to E t h = 5.5 MeV [4J. 

Radiochemical experiments integrate in time and in energy because they are slow and need time to 
produce measurable results. This causes the loss of information about single individual energy values. 
In real time experiments, instead, it is possible to obtain single values and therefore a spectrum 
energy to distinguish the different neutrino contributions. Furthermore, given that the scattered 
electron maintains the same direction of the impinging neutrino, it is possible to infer the direction of 
the incoming neutrino and therefore to point at its source. This proved that the detected neutrinos 
actually came from the Sun. The number of detected neutrinos was about 1/2 lower than the number 
of expected ones, aggravating the Solar Neutrino Problem. 

Until 1990 there were no observations of the initial reaction in the nuclear fusion chain, i.e. the 
detection of pp neutrinos, which are less model-dependent and hence more significant to test the 

lr The science that studies the interior of the Sun by looking at its vibration modes. 
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hypothesis that fusion of hydrogen powers the Sun. Two radiochemical experiments were built in 
order to detect solar pp neutrinos, both employing the reaction: 

v e + 71 G&^ 71 Ge + e- . (14.4) 

which has a threshold of E t ^ = 233 keV. 

In the Gallex experiment, located at the Gran Sasso underground laboratory in Italy, 30 tons of 
natural gallium were employed [3|6], while in the soviet-american experiment (SAGE), located in the 
Baksan underground laboratory, there were 50 tons of metallic gallium [7]. Calibration tests with an 
artificial neutrino source, 51 Cr, confirmed the efficiency of both detectors. Once again the measured 
neutrino signal was smaller than predicted by the SSM (as 60%). 



All experiments detected fewer neutrinos than expected from the SSM. Table 14.1 summarizes the 
observed vs expected ratio for all experiments. 



Homestake 


0.34 


± 


0.03 


Super-K 


0.46 


± 


0.02 


SAGE 


0.59 


± 


0.06 


Gallex and GNO 


0.58 


± 


0.05 



Table 14.1: Observed vs expected ratio in the four experiments (before SNO, see later). 



14.2.2 The advent of SNO and Kamland: the solution of the Solar Neutrino 
Problem 



The real breakthrough in solar neutrino physics was due to the advent of the SNO (Sudbury 
Neutrino Observatory) experiment. It had the peculiarity to measure simultaneously, by means of a 
deuterium Gerenkov detector, three different interaction channels for neutrinos: the neutral current 
(NC: i>x + d — > vx + p + + n), receiving contributions from all active flavors, the elastic scattering 
(ES: vx + e~ — > vx + e~) and the charged current (CC: u e + d — > e~ + p + + p + ), that is sensitive 
only to electronic neutrinos. In this way it has been possible to prove in a clear and direct way that 
the measured total neutrino flux was in very good agreement with the SSM predictions, but only a 
fraction of these neutrinos had conserved its flavor during their way from the production point in the 
Sun to the detector. 

The first SNO data jS], including elastic scattering and charged current analysis, published in 2001, 
confirmed the results obtained by previous solar neutrino experiments, mainly by Super-Kamiokande 
[9], providing a significant evidence (at the 3.3 a level) of the presence of a non-electronic active neu- 
trino component in the solar flux. For the first time it was possible to indicate the Large Mixing 
Angle (LMA) as the preferred solution of the solar neutrino puzzle, even if different alternative possi- 
bilities (and in particular the low probability, low mass -LOW- solution) were still surviving [10J. In 
the following years, the SNO experiment measured also the neutral current channel, using different 
techniques. The data of these different "phases" of the experiment are usually reported as SNO I 
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SNO II [12] (characterized by the addition of salt to improve the efficiency of neutral current 
detection) and SNO III [13J (with the use of helium chamber proportional counters). 

The year 2002 is very often denoted as the "annus mirabilis" of solar neutrino physics: on April the 
first SNO results including neutral current detection [14"] marked a turning point in the history of 
the solar neutrino problem, in October the Nobel prize for physics was awarded to R. Davis Jr. jT5j and 
M.Koshiba (for their pioneering work on the detection of cosmic neutrinos) and on December of the 
same year the first results of the Kamiokande Liquid scintillator AntiNeutrino Detector (KamLAND) 
|16j offered the first clear terrestrial confirmation of the validity of the oscillation solution to the solar 
neutrino problem. 

The total 8 B neutrino flux, <j) NC = 5.09]to;2f (stat^g^syst) x 10 6 cm" 2 s _1 , measured by SNO 
with neutral currents was in very good agreement with the SSM [17J. Assuming the standard shape for 
the component of the solar neutrino flux (undistorted spectrum hypothesis), the SNO collaboration 
recovered also a value of the non-electronic component of the flux which was 5.3 a different from zero, 
providing a direct proof of the validity of the oscillation hypothesis. These data were also decisive to 
indicate the LMA region as the solution to the solar neutrino puzzle. 

Looking at the oscillation probability it is apparent that the reactor experiments that run before 
KamLAND, and used neutrino energy beams of the order of the MeV with a baseline of the order of 
1 km, could test only values of Am 2 above 10 -3 eV 2 . The KamLAND experiment, instead, with an 
average baseline of about 180 km, was ideal to probe the LMA region, which corresponds to values of 
Am 2 of the order 10 -5 — 10 _4 eV 2 [18]. The KamLAND experiment studied the ratio of the number 
of inverse (3 decay events (due to reactor v e with an energy threshold of 3.4 MeV) to the expected 
number of events without disappearance and also the spectrum shape |16| . The observed deficit of 
events was inconsistent with the expected rate in absence of oscillation at the 99.95% confidence level. 

Since one would expect a negligible reduction of the v e flux from the SMA, LOW and vacuum solar 
neutrino solutions, the LMA was the only oscillation solution compatible with KamLAND results and 
CPT invariance. This evidence were further reinforced by the data published by the collaboration in 
the following years (with greater statistical precisions and reduced systematic errors), which showed 
also a spectral distortion in very good agreement with the oscillation solution [19|. 120] I2T] . KamLAND 
data also restricted the allowed LMA region in a significant way. The preferred values for Am 2 2 and 
#12 are slightly higher than the ones corresponding to the best fit solution of the solar neutrino exper- 
iments, but this small tension can be explained by taking into account the experimental uncertainties. 
Moreover, the difference on the Amf 2 parameter has been reduced by the more recent solar neutrino 
data. 



14.3 Standard Solar Model 

SSMs have to be understood, primarily, as a framework within which solar models can be constructed 
and clear predictions can be made with respect to the properties of the solar interior, including the 
production of solar neutrinos. The defining characteristics are simple: the SSM is the result of the 
evolution of a 1 M star since its formation and, the evolutionary models have to include the physical 
ingredients considered standard in stellar structure and evolution models (here, standard also implies 

2 For instance, in a simple 2 flavor analysis, the flavor transition probability is given by the expression 
P\2 = sin 2 (2612) sin 2 ( — - J > where #12 is the mixing angle between the two flavors, Amjj = m\ — m 2 , the 
difference of the masses squared, L the distance traveled, and E the neutrino energy. 
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trying to keep to a minimum the number of free tunable parameters -knobs- in the model). SSMs are 
therefore progressively refined as our understanding of stellar physics progresses. 

In practice, a SSM is constructed as follows. An initial chemically homogeneous model of a 
1 Mq stellar model on the pre-main sequence is constructed with a composition determined by a 
guess (educated one) for the initial mass fractions of hydrogen X- m \, helium Y m \, and metals Z- m \ 
{X- m \ + lini + Z- m \ = 1); additionally, a third free parameter has to be specified, the mixing length 
parameter ckmlt of convection. This model is then evolved up to the solar system age r = 4.57 Gyr 
|22t [23] . At this age the model is required to match the present-day solar luminosity L and radius 
Rq, as well as the surface metal-to-hydrogen abundance ratio (Z/X) Q . The initial and final surface 
metal-to-hydrogen ratios differ by about 10 to 15% due to the effects of gravitational settling. In 
general, the SSM constructed with the first set of guesses for omlt, ^Ini, and Z- m \ will not lead to a 
satisfactory agreement with the surface constraints, and an iterative procedure is used to refine the 
free parameters until the right surface conditions are achieved at tq. In general, surface conditions are 
matched to one part in 10 5 or 10 6 within two or three iterations. It is important to keep in mind that 
the SSM is not just a snapshot aimed at representing the present-day structure of the Sun, but actually 
the result of taking into account all its previous history. There are alternative ways to construct a 
model of the present-day solar structure using, for example, helioseismic constraints. These kind of 
models are constructed 'ad-hoc' to match helioseismic data and have, therefore, limited predictive 
power. 

The internal structure of a SSM depends on the values adopted for the three constraints mentioned 
above and, of course, on the physical inputs of the models such as the radiative opacities, cross sections 
of nuclear reactions and others. Next we describe the changes/updates that have occurred during the 
last decade that impact predictions of solar models. 

14.3.1 Input physics and parameters 
Solar Surface Composition 

The constraint imposed by the surface metallicity of the Sun or, more precisely, the surface metal-to- 
hydrogen ratio (Z/X) & , is critical in the construction of solar models. The reason is that, aside from 
the 10 to 15% change in this value due to the action of gravitational settling, (Z/X) & determines 
almost directly the metallicity of solar models. As for any other star, the metal content in the Sun 
has a fundamental role in its structure through its contribution to the radiative opacity k, which 
determines, in turn, the temperature gradient in the radiative solar interior. It is important, in 
fact, that the abundance of individual metals are accurately determined, because different elements 
contribute to the radiative opacities in different regions of the Sun. 

The abundance of metals in the solar surface has to be determined or inferred from a variety 
of sources: photospheric abundances from solar spectra, chemical analysis of primitive meteorites, 
emission lines from the solar corona, composition of the solar wind [53] . While meteoritic abundances 
are the most precisely determined, at least 2/3 of the solar metallicity is composed by the volatile 
elements C, N, and O and can only be determined from analysis of the solar spectrum. 

Over the last decade, the development of three-dimensional radiation hydrodynamic (3D RHD) 
models of the solar atmosphere has prompted a thorough revision of the solar composition determined 
from the solar spectrum. These 3D RHD models of the solar atmosphere capture the dynamics of 
convection and its interaction with the radiation field, and are able to reproduce features such as the 
solar granulation pattern, observed limb-darkening, asymmetries in the shapes of spectral lines [25]. 
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The structure of the solar model atmospheres derived by different groups are nicely consistent with 
each other, adding to the credibility of the models. Newly derived spectroscopic abundances rely on 
the 3D atmosphere model, or more appropriately on a one- dimensional model obtained from a suitably 
averaged 3D model, as the background on top of which detailed radiative transfer and line formation 
calculations are performed a posteriori. It is this second step that leads, finally, to the determination 
of the abundances of the different elements. The most thorough and consistent determination of the 
solar photospheric abundances based on 3D model atmospheres has been presented by Asplund and 
collaborators [261 [27] , although revision on key elements like oxygen were initially published already 
in 2001 [28j. In addition to using 3D RHD atmosphere models, non-local thermodynamic equilibrium 
has been taken into account when computing line formation for some key elements such as C, N, and 
O. Also, and this is of particular importance for oxygen, blends in the solar spectrum that had been 
previously unnoticed were identified and taken into account in the determination of abundances. The 
most relevant result in the context of solar models and neutrinos is that abundances of CNO elements 
(also Ne, but this is mostly because its abundance ratio to oxygen is assumed fixed) have been revised 
down by 30 to 40 %. Combining the abundance of all metals, the present-day metal-to-hydrogen ratio 
that has been obtained is (Z/X) & = 0.0178 |26j . This represents a large decrease in comparison with 
previously accepted values, 0.0245 [29J and 0.0229 [30], that have been widely used in solar modeling. 
We note, however, that results by Asplund have not been unchallenged. In fact, also based on 3D RHD 
model solar atmospheres, larger CNO abundances have been derived [3l] to yield (Z/X) Q = 0.0209, 
much closer to older determinations. Discrepancies between authors seem to have their origin at the 
preferred set of spectral lines each group uses and on using either a spectral synthesis or equivalent 
width techniques to determine the final abundances. 

In the last decade there have been two flavors in SSM calculations. In one case a high solar 
metallicity from older determinations |29| [30] is adopted; we will generically refer to these models as 
high-Z solar models. In the other case a low (Z/X) & |26t [2"7] is taken from and we refer to these, not 
surprisingly, as the low-Z solar models. Differences in the structure of high-Z and low-Z models are 
readily noticeable in quantities such as the internal sound speed and density profiles, the depth of the 
solar convective envelope, and the surface helium abundance among others. The deficit that low-Z 
models have in matching helioseismic constraints has been named the solar abundance problem in the 



literature, in clear analogy to the solar neutrino problem. We discuss it in some detail in Section 14.3.2 



Radiative Opacities 

The most widely used calculations of atomic radiative opacities, appropriate for solar interiors, are 
those from OPAL [32J. However, the Opacity Project (OP) released in 2005 a completely independent 
set of atomic radiative opacities for stellar interiors [33J. In the case of the solar radiative interior, 
differences between OPAL and OP Rosseland mean opacities are of the order of a few percent, with 
OP being larger by about 3% at the base of the convective zone and 1 to 2% smaller in the central 
regions (see Fig. 7 in [33] ). At low temperatures, at which molecules can form, neither OP or OPAL 
atomic opacities are adequate and have to be complemented by low-temperature opacities |34| . Due 
to the relatively high solar temperature, their influence in the properties of solar models is rather 
limited. 
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Nuclear reactions cross sections 

Experimental and theoretical work on the determination of nuclear cross sections have been very 
active fields with a strong impact on solar model predictions of solar neutrino fluxes. Recently, a set 
of recommended rates and uncertainties, expressed through the S-factoiJ^J for all the reactions both 
in the pp-chains and CNO-bicycle that are relevant to solar modeling and neutrino production, has 
been published (Solar Fusion II, [35], hereafter SFII). The results presented in SFII reflect the progress 
made in laboratory and theoretical nuclear astrophysics over the last decade, since the publication 
of the seminal Solar Fusion I (SFI) article [36] . Unfortunately, for reasons of space, here we cannot 



review in detail every reaction. Instead, we provide in Table 14.2| the standard S'-factors at zero energy, 
S(0) , and the uncertainties recommended in SFII for the most relevant reactions. For comparison, 
with results from SFI are also shown. The impact of changes in key reactions on the production of 



neutrino fluxes is discussed in Section 14.3.3. The reader is referred to the SFII paper and references 
therein for details on the experimental and theoretical developments in nuclear astrophysics related 
to the Sun during the last decade. 

Reaction SFII SFI 

5(0) [keV b] 5(0) [keV b] 

~Su P{P, e + v e )d 4.01 x 1(T 22 (1± 0.010) 4.00 x 10~ 22 (1 ±0.005) 

5 33 3 He( 3 He, 2p) 4 He 5.21 x 10 3 (1 ± 0.052) 5.4 x 10 3 (1 ± 0.074) 

5 34 3 He( 4 He, 7) 7 Be 5.6 x lO" 1 (1 ± 0.054) 5.3 x lO" 1 (1 ± 0.094) 
S hep 3 He(p, e+^ e ) 4 He 8.6 x 10" 20 (1 ± 0.30) 2.3 x 10~ 20 

Sit 7 Be(p, 7 ) 8 B 2.08 x 10~ 2 (1 ± 0.077) 1.9 x 10~ 2 (l±oio) 
Si,i4 14 N(p, 7 ) 15 1.66(1 ±0.072) 3.5 (l±g;^) ' 

Table 14.2: Standard astrophysical factors and uncertainties for key nuclear reactions in the pp- 
chains and CNO-bicycle. SFII represents the state-of-the-art [35]; SFI [36] shows, for comparison, the 
situation around 1998. 



14.3.2 Solar Models: Helioseismology 

Helioseismology, the study of the natural oscillations of the Sun, provides a unique tool to determine the 
structure of the solar interior. The '90s witnessed a rapid development of helioseismic observations and 
analysis techniques, which led, in very few years, to an accurate characterization of the solar interior 
[37] • The agreement between SSMs and helioseismic inferences of the solar structure |38| I17| provided 
a strong support to the accuracy with which SSMs could predict the 8 B neutrino flux and, therefore, 
a strong indication, before Kamland and SNO results found evidence of neutrino flavor oscillations, 
that the solution to the solar neutrino problem had to be found in the realm of particle physics. 

In the context of the present article, the most relevant results from helioseismology are the follow- 

3 A non-resonant charged-particle induced reaction cross section can be written as cr[E) = S ^ exp [— 2-7r77(E)] where 
rj(E) = ZiZ^a/v is the Sommerfeld parameters, v — ^/\2E/[i), a the fine structure constant in natural units, and [i the 
reduced mass of the interacting nuclei. The nuclear physics is isolated in S(E), the astrophysical or S-factor, a slowly 
varying function of energy that can be more accurately extrapolated from experimental data down to the energy of the 
Gamow peak. 
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ing. The depth of the convective envelop^] is Rqz = 0.713 ± 0.001 R© [39j and the surface helium 
abundance Y$ = 0.2485 ± 0.0034 [10]. The sound speed differences between the Sun and a reference 
solar model can be obtained by inversion from the oscillation frequencies with a formal error of a few 
parts per 10~ 4 for most of the solar interior 0.07 < R/Rq ^ 0.95 [JIJ H2]. Most recently, using a 
time series 4752 days-long from the Birmingham Solar Oscillation Network, improved results on the 
sound speed in the solar core have been obtained [IB]. The density profile can also be determined from 
inversion of frequencies, but with worse precision than for the sound speed, and we therefore assign 
to it a secondary role in constraining the solar structure. 

As mentioned previously, metals determine to a large extent the radiative opacity in the solar 
interior and, in this way, define the temperature stratification from below of convective envelope 
inwards, to the solar center. At the base of the convective zone, for example, metals are responsible of 
about 70% of the total radiative opacity with O, Fe and Ne being the main contributors. In the solar 
core, where light metals are completely ionized, the contribution from Fe and, to a lesser extent Ni, 
Si and S, is still above 30%. In view of this, it is not surprising that the low CNO and Ne abundances 
determined from 3D model atmospheres have a strong impact in the structure of the solar interior. 

It has been clear since initial works where low-Z SSMs were presented that low (Z/X) & values 
posed a problem, later named the solar abundance problem, for solar modeling \4A\ 1131 HOl 06] . In 
short, all helioseismic predictions of these models are in disagreement with observations. On the other 
hand, high-Z SSMs have consistently reproduced earlier success [38 . The solar abundance problem 
represents the incompatibility between the best solar atmosphere and interior models available }47j . 
In this review, we will base the presentation and discussion of results on the most up-to-date standard 
solar models that we identify as SFII-GS98 and SFII-AGSS09 |48j . representative of high-Z and low-Z 



SSM families defined in Section 14.3.1 respectively. With the exception made on small quantitative 
variations, results based on these models are extensible to results for all SSMs available in the literature 
corresponding to each of the two families. 

The most important characteristics of the SFII-GS98 and SFII-AGSS09 models are summarized 



in Table 14.3 Helioseismic constraints are also included for comparison when appropriate. The 
disagreement between SFII-AGSS09 and helioseismic data is evident in the surface metallicity and 
helium abundances, Z$ and Yg, and in the depth of the convective envelope Rcz- When model 
uncertainties are included, the discrepancy between SFII-AGSS09 and seismic results are, for each of 
the quantities mentioned above, of the order 3 to 4 — a [49J. On the contrary, the SFII-GS98 model 
performs very well, within 1 — a when model uncertainties are accounted for. 



Very explicit manifestations of the solar abundance problem are shown in the plots in Figure 14.1 
where degradation in the sound speed and density profiles found in low-Z SSMs are clearly evident. 
Particularly the peak in the sound speed profile differences found right below the convective zone is 
4 times larger in the low-Z SFII-AGSS09 than in the high-Z SFII-GS98 model. The reason is the 
wrong location of Rcz m the model, caused by the lower opacity which, in turn, is due to the low 
abundance of metals. The density profile also shows very large discrepancies, but they are less telling. 
Density inversions include as a constraint the known value of the solar mass and for this reason small 
differences in the core, where density is large, translate into the large difference seen in the outer 
envelope. The average rms in the sound speed and density differences, (5c/ c) and (5p/p), also show 



4 The Sun is characterized by an outer region where energy is transported by convection. The boundary between this 
region, located at Rcz, and the radiative interior can be accurately located by helioseismology because the discontinuity 
in the slope of the temperature gradient accross this boundary leaves its imprint in the solar sound speed profile. The 
depth of the envelope can be located by helisoeismology because properties of solar oscillations are sensitive to the 
derivative of the sound speed as a function of depth. 
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SFII-GS98 


SFII-AGSS09 


Helioseismology 


(Z/X) Q 


0.0229 


0.0178 






0.0170 


0.0134 


0.0172 ± 0.002 |50J 


Y S 


0.2429 


0.2319 


0.2485 ± 0.0034 @U] 


-Rcz/R© 


0.7124 


0.7231 


0.713 ±0.001 [39] 


(Sc/c) 
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Table 14.3: Main characteristics of SSMs representative of high-Z (GS98) and low-Z (AGSS09) solar 
compositions. Models have been computed including the most up-to-date input physics [48]. Helio- 
seismic constraints are given when available. See text for details. 



that low-Z models are about 4 times worse than high-Z models. 
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Figure 14.1: Sound speed and density relative differences between solar models and the Sun as de- 
termined from helioseismic inversions |49j . The convective envelope is depicted by the grey area. 



Low-degree helioseismology provides useful information about the solar innermost regions. Specific 
combinations of mode frequencies enhance the signal that the structure of the solar core imprints on 
the oscillation pattern |52j . This has been used to determine the mean molecular weight averaged 
over the innermost 20% solar core [HI], (mc) m Table 14.3 Comparison with SSMs results shows that 
{fj-c) is too low in low-Z models as a result of the lower helium abundance (Yc)- This is due to the 
lower temperature in the solar core and the constraint imposed by the solar luminosity. The decreased 
nuclear energy production originated by a smaller core temperature has to be compensated by an 
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increased hydrogen mass fraction, therefore leading to a lower molecular weight. It is interesting to 
note this puts a stringent constraint in the amount of rotational mixing that can take place in the 
solar core if the low-Z abundances are correct, since any mixing would lower the molecular weight even 
more, by bringing fresh hydrogen from outer regions, and make the agreement with helioseismic data 
worse. 

The current situation regarding SSMs and their performance against helioseismic inferences on 
the solar structure can be summarized as follows. SSMs that use solar abundances derived from 
ID model atmospheres [29, 30J, i.e. high-Z models, reproduce overall the most important seismic 
constraints. Improvements in the input physics, e.g. radiative opacities and nuclear reaction rates, 
that have occurred over the last 10 years introduce only small changes to the solar structure as seen by 
helioseismology. On the other hand, the solar abundance problem arises if the solar surface composition 
used to construct SSMs are derived from the most sophisticated 3D RHD solar model atmospheres. 
The family of low-Z SSMs does not match any helioseismic constraint. 

Have we reached the limit where the paradigm of the SSM is not good enough as a model of the 
solar interior? Are the 3D-based determinations of solar abundances systematically underestimating 
the metallicity of the solar surface? Does the microscopic input physics in solar models, e.g. radiative 
opacities, need to be thoroughly revised? It is not possible to advance answers to these questions, but 
solar neutrino experiments can play an important role in guiding research towards the solution of the 
solar abundance problem. In the next section we discuss the current status on the theoretical predic- 
tions of solar neutrino fluxes and the prospects of using solar neutrinos to constraint the properties of 
the solar core. 



14.3.3 Solar Models: Neutrino Fluxes 
Production 

Based on theoretical arguments and indirect evidence, it has long been believed that the source of 
energy of the Sun is the conversion of protons into helium, 4p — > 4 He + 2e + + 2u e + 7. The original 
quest for solar neutrinos was indeed the search for the experimental confirmation of this hypothesis. 
In more detail, hydrogen burning in the Sun (and in all other hydrogen-burning stars) takes place 
either through the pp-chains or the CNO-bicycl^] |53| [23] . Proton fusion through the pp-chains is a 
primary process because only protons need be present in the star. On the contrary, the CNO-bicyle 
is secondary because proton fusion relies on, and is regulated by, the abundance of C, N, and O 
which act as catalyzers. This qualitative difference is very important, since it renders neutrino fluxes 
from the CNO-bicyle a very good diagnostic tool to study properties of the solar core, particularly its 
composition, as it will be discussed below. A general discussion on the production of solar neutrinos 
is out of the scope of the present review, but can be found elsewhere [23J . 

SSM calculations of neutrino fluxes have been affected by developments in the input physics dis- 
cussed in previous sections. The two areas that have the strongest impact on the neutrino fluxes 
predicted by models are: changes in nuclear cross sections, and the new solar composition. In Ta- 



ble 14.4 we list the results for neutrino fluxes for the up-to-date SSMs SFII-GS98 and SFII-AGSS09. 



For comparison we include, in the last column, results from the BP04 SSM |54j . 

Under peculiar conditions reached in advanced phases of stellar evolution, hydrogen can be converted into helium 
by other cycles like the NaMg-cycle. While important for nucleosynthesis or intermediate mass elements, these processes 
are not energetically relevant 
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pep 

hep 

7 Be 

8 B 


5.98(1 ±0.006) 
1.44(1 ±0.012) 
8.04(1 ±0.30) 
5.00(1 ±0.07) 
5.58(1 ±0.13) 


6.03(1 ±0.006) 
1.47(1 ±0.012) 
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4.56(1 ±0.07) 
4.59(1 ±0.13) 
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7.8(1 ±0.16) 
4.86(1 ±0.12) 
5.79(1 ±0.23) 


13 N 


2.96(1 ±0.15) 


2.17(1 ±0.13) 


< 6.7 


5.71(1 ±0.36) 


15Q 


2.23(1 ±0.16) 


1.56(1 ±0.15) 


< 3.2 


5.03(1 ±0.41) 


17 F 


5.52(1 ±0.18) 


3.40(1 ±0.16) 


< 59. 


5.91(1 ±0.44) 


x 2 /P agT 


3.5/90% 


3.4/90% 







Table 14.4: SSM predictions for solar neutrino fluxes (second and third columns) and solar neutrino 
fluxes (fourth column) inferred from all available neutrino data. Units are, in cm _2 s _1 , as usual: 10 10 
(pp), 10 9 ( 7 Be), 10 8 (pep, 13 N, 15 0) 10 6 ( 8 B, 17 F), and 10 3 (hep). 
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Figure 14.2: Normalized production profiles of solar neutrinos as a function of solar radius. 



The most striking difference is the large reduction in the 13 N and 15 fluxes between the SFII- 
GS98 and BP04 models, which use the same solar composition. This reduction comes as a result of 
the new determination of Si 14, mostly by the LUNA experiment [551 156] , that has halved its value 
with respect to previous results (Table 14.2). If correct, the new expectation value of the combined 
13 N± 15 fluxes poses an even more challenging task for neutrino experiments to detect CNO fluxes. 
By comparing fluxes in Table 14.4 for models computed with the same solar composition (SFII-GS98 
and BP04), it can be seen that in terms of flux values, those associated with the pp-chains have not 
changed much since 2004, despite improvements in the input physics entering solar model calculations. 
Few percent changes are present and are the result of changes in the nuclear cross sections discussed 
before and also of the new OP radiative opacities. This is an encouraging situation; it implies that 
neutrino fluxes are robust predictions of solar models and, as experimental data on solar neutrinos 
accumulate, it will be possible to start fulfilling the initial goal posed by Davis and Bahcall: to use 
solar neutrinos to learn about the solar interior. 



In Figure 14.2 we show the distribution of the solar neutrino fluxes as a function of solar radius. 
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Together with the electron density profile, provided also by solar models (and neutron density profiles 
for sterile neutrino studies), these quantities are of fundamental importance for neutrino oscillation 
studies. It is worth noting that the 13 N flux has two components. The larger one is associated with 
the operation in quasi steady-state of the CN-cycle in the innermost solar core (R < 0.1 R©), and for 



this reason coincides with the production region of the 15 flux (right panel in Figure 14.2, blue and 
black curves respectively). This component of the 13 N flux, as well as the total 15 flux, is linearly 
dependent on Si 14. The additional component of the 13 N flux comes from the residual burning of 12 C 
by the reactions 12 C(p, 7) 13 N(/3 + ) 13 C at temperatures not high enough to close the CN-cycle with a 
proton capture on 14 N. This component is completely independent of Si. 14. The careful reader will 
notice that the ratio of 13 N and 15 fluxes is different in the SFII-GS98 and BP04 models, despite 
having the same solar composition. Whereas the added 13 N+ 15 neutrino flux is linearly proportional 
to the C+N abundance in the solar core and also linearly proportional to 81,14, this degeneracy can 
be broken, at least theoretically, if the two fluxes can be experimentally isolated from one another. 
The impact of the low-Z solar composition on the production of solar neutrinos can be grasped 



by comparing results of models SFII-GS98 and SFII-AGSS09 shown in Table |14.4| As stated before, 
metals shape the solar structure through the radiative opacity. The lower abundance of metals in the 
AGSS09 composition is responsible for a reduction of the temperature in the solar core of about 1%. 
Because of the extreme temperature sensitivity of some of the neutrino fluxes this is enough to produce 
large changes in the total fluxes. The most extreme case is, of course, 8 B, with the SFII-AGSS09 value 
being ~ 20% smaller. For 7 Be the reduction is of ~ 9%. 

Given the small uncertainties in the experimental determination of these fluxes, it would be tempt- 
ing to think these neutrino fluxes have the potential to discriminate between the two flavors of solar 
composition and contribute, in this way, to the solution of the solar abundance problem. As can be 



seen in Table 14.4, unfortunately, the 7 Be and 8 B fluxes determined from experiments lie almost right 
in between the high-Z and low-Z models. 

In any case, since it is known that low-Z solar models do not reproduce well the solar structure 
as discussed in the previous section, it is dangerous to extract conclusions from comparing neutrino 
fluxes of this model to experimental results. Regardless of what the solution to the solar abundance 
problem is, since it will modify the solar interior structure, it will also change the expected values 
for the neutrino fluxes. In this regard, CNO fluxes are particularly interesting. Although they are of 
course affected by temperature variations to a comparable degree as the 8 B flux is, they carry an extra 
linear dependence on the solar composition that is not related to temperature variations. Of particular 
interest is the linear dependence of the 13 N and 15 fluxes on the combined C+N abundance^] It is 
this dependence that enhances their capability as a diagnostic tool. In fact, differences between SFII- 
GS98 and SFII-AGSS09 models for these two fluxes are of the order of 30% (taking SFII-GS98 as 
reference) and, what is more important, a large contribution to these differences does not have an 
origin on temperature differences between the models. 

The last row in Table 14.4 shows the results of a x 2 test f° r the two models against the solar 



fluxes also shown in the table. It is clear that both SSMs give very good agreement with current data. 
We emphasize again, however, that the four fluxes that are currently well determined from data and 
the luminosity constraint, depend on the solar composition only in an indirect manner. Experimental 
determination of the combined 13 N+ 15 flux will therefore provide qualitatively new information 
on the solar structure and composition. In fact, one can take advantage of the similar response to 



6 The 17 F flux is linearly dependent on O, but unfortunately the flux is too low to be detectable with current experi- 
mental capabilities. 
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temperature variations that CNO fluxes and the 8 B flux has. This has been exploited [57] to develop 
a very simple method to determine the solar core C+N abundance that minimizes environmental 
uncertainties in solar models (that is, sources of uncertainty that affect the solar core temperature). 
The idea is simple: the temperature dependences are cancelled out by using an appropriate ratio 
between the 8 B and the combined 13 N+ 15 fluxes where SSM fluxes only act as normalization values 
and the overall scale is determined by an actual 8 B flux measurement. The only additional requirement 
is that a measurement of the combined 13 N+ 15 flux becomes available. The current upper limit on 
this combined flux from Borexino |58j places an upper limit on the C+N central mass fraction of 
X c+ n < 0.072. Results for the SFII-GS98 and SFII-AGSS09 models are X c+N = 0.048 and 0.039 
respectively. 



Uncertainties 



Uncertainties in the model predictions of solar neutrino fluxes are given in Table 14.4 For deriving the 
total uncertainty basically two approaches can be used. On one hand, all contributions of uncertainty 
can be treated simultaneously by doing a Monte Carlo simulation [59J. The advantage is that intrinsic 
non-linearities are captured in the total error. The disadvantage is that individual contributions to 
the total uncertainty are hardwired in the final result and can not be disentangled. Fortunately, 
for the current level of uncertainties entering SSM calculations, non-linearities seem to be negligible 
and the total uncertainty in neutrino fluxes can be obtained (adding quadratically) from individual 
contributions. To compute the latter, the expansion of fluxes as a product of power-laws in the input 
parameters [23j around central values is a widely used, practical, insightful and accurate approach. 



Uncertainties in the model fluxes listed in Table 14.4 have been obtained in this way. 



The most important change introduced in the estimation of uncertainties is related to the treatment 
of the solar composition. Up until the BP04 model |54| . the uncertainty in the solar composition was 
taken into account by considering variations of the total solar metallicity (to be more precise, changes 
in the (Z/X) Q value used to construct SSMs). This leads to an overestimation of the neutrino 
uncertainties. The reason is that metals dominating the error budget in {Z/X) & (C, N, O, and 
Ne) have, at most, a moderate impact on the neutrino fluxes because of their small contribution 
to the radiative opacity, and therefore a rather small impact on temperature, in the region where 
most neutrinos are produced. On the other hand, elements such as Fe, S and Si are second order in 
determining (Z/X) & but play a fundamental role as sources of opacity in the solar core. It is important, 
therefore, to treat metal uncertainties individually [60]. Of course, in the case of the CNO fluxes the 
situation is different because CNO elements catalyze the CNO-bicycle and this overimposes an almost 
linear dependence of the 13 N and 15 on the C+N content of the solar and a similar dependence 
of 17 F on the O abundance. The uncertainties in the neutrino fluxes given for the SFII-GS98 and 
SFII-AGSS09 SSMs have been computed using the uncertainties for each relevant element given in 
the original publications [60\ 126]. As a result, for either family of solar models, i.e. high- Z or low-Z 
models, the solar composition is not the dominant source of uncertainty for any of the fluxes of the 
pp-chains. In the case of the CNO fluxes, the linear dependence mentioned above is the dominant 
source of uncertainty: the combined C+N abundance contributes to a 12% uncertainty for both the 
13 N and the 15 fluxes, and the O abundance to 15% in the 17 F flux. 

In the case of the non-composition uncertainties, the situation has improved in some cases thanks 
to more precise measurements of nuclear reaction rates. This is the case, in particular, for the 
3 He( 4 He, 7) 7 Be reaction, which now contributes only 4.7% and 4.5% of the total uncertainty in the 
7 Be and 8 B fluxes respectively. For comparison, the analogous contributions in the BP04 model were 
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Table 14.5: Percentage contribution of selected individual sources of uncertainty to the neutrino fluxes. 



8.0% and 7.5% Significant progress has also been achieved regarding 14 N(p, 7) 15 0, which now 

introduces uncertainties of only 5% and 7% in the 13 N and 15 O fluxes, half the amount it did in 2004. 
An important contribution to the uncertainty in the 8 B flux now comes from 7 Be(p, 7) 8 B because the 
uncertainty of this reaction has been revised upwards [35] . Even if the uncertainty in this rate is now 
smaller than in SFI (see Table 14.2), it is larger than that used for the BP04 model, which was taken 
considering only one experimental result for this reaction. 

While progress has been done in some cases, others have not seen much development, particularly 
diffusion and the delicate issue of radiative opacities. In Table 14.5 we give the individual contributions 
to flux uncertainties for the most relevant sources. The reader can compare directly to the situation 
in 2004 [Ml- 



14.4 Neutrino flavor conversion in vacuum and matter 

Neutrino flavor conversion has been reviewed by A. Yu. Smirnov in this volume and we refer the 
reader for a detailed physics discussion and references to his article. Here we just summarize the basic 
features and formulae of flavor conversion relevant to solar neutrinos. 

We consider mixing of the three flavor neutrinos. The description of flavor conversion of solar 
neutrinos traveling through a medium is simplified because a) the hierarchy in mass splittings deter- 
mined by solar and atmospheric data leads to a reduction of the three neutrino flavor conversion to 
an effective two flavor problem, b) the neutrino parameters, the mixings and solar mass splitting, lead 
to adiabatic flavor conversion in solar matter and to cancel the interference term by averaging out. 
Therefore, the physics of the flavor conversion of solar neutrinos is described by simple expressions 
with very good accuracy. In practice, the survival probability is computed numerically to correctly 
include the number density of scatterers along the trajectory of neutrinos from production to detection 
and to average over the neutrino production region. In solar neutrino flavor conversion, Uu and v T 
are indistinguishable and therefore the survival probability of electron neutrinos is the only function 
needed to describe the flavor composition of the solar neutrino flux. 

Solar neutrino survival or appearance probabilities depend on three oscillation parameters: the 
solar oscillation parameters (#12, Ato^), and #13. The survival probability in the absence of Earth- 
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matter effects, i.e., during the day, is well described by 



p Fe = cos" y 13 



- + - • 005 26*5 • 008 26*12 j -!- sin f) l: >. 



Here 9s is the mixing angle at the production point inside the Sun: 

cos 26*5 = cos26 m (ps) 
where 9 m (p) is the mixing angle in matter of density ps, 

cos 2#12 — £s 



cos 26*5 
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(14.5) 
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In (|14.7|), £s is defined as the ratio of the neutrino oscillation length in vacuum, l v , to the refraction 
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In (14.8) and (14.9), ps is the solar matter density, Y e $ is the number of electrons per nucleon, and 
is the nucleon mass. The electron solar density and neutrino production distribution of the neutrino 
fluxes are derived from solar models as discussed in previous section. In the last line in (14.8) we have 
used the best fit values of the global analysis Am 2 = 7.5 x 10 -5 eV 2 . The ratio of the parameter ps to 
008 26*12, separates the region where the flavor conversion corresponds to vacuum averaged oscillations 
from the one of matter dominated conversion. 

The v e survival probability at night during which solar neutrinos pass through the Earth can be 
written as 



P ee = P ee ~ cos 29 S COS 2 6>i 3 (/ re9 ) ze nith 



(14.10) 



where is the one given in (14.5). f reg denotes the regeneration effect in the Earth, and is given as 
freg = P 2e — sin 2 6*i2 cos 2 #13, where P2e is the transition probability of second mass eigenstate to v e . 



Under the constant density approximation in the Earth, f reg is given by 



freg = £e cos 2 6*13 sin 2 29 E sin 2 



a E COS 2 0i 3 (l - 2C E l COS 2 012 + ££ 2 )5 



L 



(14.11) 



for passage of distance L, where we have introduced a E = V2G E nf arth = y2^£mXse _ 



In (14.11), 9 E and £ E stand for the mixing angle and the ^ parameter [see (14.8)] with matter 
density p E in the Earth. Within the range of neutrino parameters allowed by the solar neutrino data, 
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the oscillatory term averages to | in a good approximation when integrated over zenith angle. Then, 
the equation simplifies to 

(/reg) zenith = ^ COS 2 Ol^E SU1 2 28 E . (14.12) 

At E = 7 MeV, which is a typical energy for 8 B neutrinos, £ E = 3.98 x 1(T 2 and sm26 E = 0.940 
for the average density p E = 5.6g/cm 3 and the electron fraction Y eE = 0.5 in the Earth. Then, 
(/reg) zenith is given as (/reg) z enith = 1-76 x 10~ 2 for the best fit neutrino parameters. This result 
is in reasonable agreement with the computed Earth-matter factor using the best estimates on the 
Earth-matter density. 

14.5 Recent solar neutrino measurements 
14.5.1 The SNO and SK legacy 

After the results and analyses from 2002, it was clear that the LMA oscillation was the right 
solution of the long standing solar neutrino puzzle [61 J , but the activity of the SNO and SK experiments 
continued in the following years. The data obtained from these experiments were very important in 
making the LMA solution more robust and in improving the accuracy and precision of the mixing 
parameters determination. 

The so-called SNO II experiment began in June of 2001 with the addition of 2000 kg of NaCl 
to the 1000 tons of D2O and ended in October 2003 when the NaCl was removed. The addition of 
salt significantly increased SNO's efficiency (by a factor ~ 3 with respect to the pure D2O phase) in 
the detection of neutrons produced in the neutral current (NC) disintegration of deuterons by solar 
neutrinos and, by enhancing the energy of the 7-ray coming from neutron capture, allowed a more 
precise measurement of this interaction channel, well above the low-energy radioactive background. 
Moreover, the isotropy of the multiple 7-ray emission by neutron capture on 35 CI is different from the 
one of the Cerenkov light emitted by the single electron of the charged current interaction; therefore, 
by studying the event isotropy, it has been possible to separate the neutral from the charged current 
events without any additional assumption on the neutrino energy spectrum. The salt phase results 
have been reported in two main publications. In |12j . referring to the first 254 live days, a global 
analysis including all the solar and reactor neutrino results rejected the maximal mixing hypothesis at 
a 5.4 a level and gave a value of the 8 B neutrino flux in agreement with previous measurements and 
with SSMs. These results were essentially confirmed (even if with a small shift towards larger values 
of the mixing angle) by the second publication [62], which included the full data of the salt phase (391 
live days), analyzed in terms of the CC spectra (starting from 5.5 MeV kinetic energy) and NC and 
ES integrated fluxes separately for day and night. The day- night asymmetry in the neutral current 
rate, which would be an indication of oscillation to sterile neutrinos or non standard interaction with 
matter in the earth, came out to be consistent with zero. 

This result confirmed also the outcome of the study performed for elastic scattering (ES) inter- 
action above 5 MeV by the Super-Kamiokande collaboration |63] . The full SK-I low energy data, 
corresponding to 1496 live days until July 2001, were investigated analyzing the time variations of 
the ES rates and fitting them to the variations expected from active two neutrino oscillations. The 
day-night asymmetry turned out to be Adn = ^pinp = —0.021 ± 0.020 (stat.) lo!oil ( s y s *-)> which 
is consistent with zero within 0.9 a. This value was in good agreement also with the LMA oscillation 
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solution, which (for the best fit parameter) predicted [63J Aon = —0.018 ± 0.016 (stat.) +o 012 ( s y s t-) 



-0.012 

The SK analysis [63^ [64] also showed that the energy spectrum of the recoiling electron was consis- 
tent with an undistorted solar 8 B neutrino spectrum and did not find any anomalous periodic time 
variation of the rates, apart from the expected seasonal variation due to the Earth's orbit eccentric- 
ity. The SK best fit point was in quite a good agreement with the SNO results, even if SK would 
favor slightly larger values of tan 2 #. A SNO-only analysis gave the following best fit parameters [62]: 
Amj 2 = 5.0 x 10~ 5 eV 2 , tan 2 #i 2 = 0.45. Including all the other solar neutrino and the KamLAND 
results the best fit was obtained for Amf 2 = 8.0l^ x 10~ 5 eV 2 , tan 2 6>i 2 = 0.452^°;°^. The effect of 
KamLAND data was mainly to increase the value of Am 2 and to restrict the allowed region in the 
mixing parameter plane. The main difference of the global analysis done with the SNO salt phase 
data with respect to previous studies was the possibility to exclude at 95%C.L. the secondary region 
at even larger values of the mass differences (the so called LMA II solution, with Am 2 2 > 10 _4 eV 2 ). 

In the third SNO phase (November 2004-November 2006) the neutral current signal neutrons 
were mainly detected by means of an array of s He proportional counters deployed in the Z) 2 and 
looking at the gas ionization induced by neutron capture on 3 He. In this way the fluxes correlation 
was reduced and the accuracy in the mixing angle determination was improved. The total active 8 B 
neutrino flux was found [13] to be 5.54+Q3i(stat)+Q34(syst) x 10 6 cm~ 2 s -1 , in agreement with previous 
measurements and SSMs. The ratio of the 8 B neutrino flux measured with CC and NC reaction was 
^cc°/^nc° = 0.301 ±0.033. The global solar neutrino experiment analysis included, in this case, also 



the first results coming from the Borexino experiment [65], that we discuss in subsection 14.5.3 The 
best fit point moved to Am 2 2 = 4.90 x 10 _5 eV 2 , tan 2 #i 2 = 0.437 and the uncertainty in the mixing 
parameter plane was still quite large. Adding the KamLAND data, the allowed region was significantly 
restricted (mainly for Am 2 ) and the marginalized la regions were Amf 2 = 7.59+ 

oil x 10 5 eV 2 , 

tan 2 12 = 0.469ig;g£. 

A subsequent joint reanalysis of SNO I and SNO II data, known as LETA (Low Energy Threshold 
Analysis) |66j . succeeded, with improved calibration and analysis techniques, in lowering the energy 
threshold, with respect to previous analyses ([671 E2]), down to an effective electron kinetic energy 
of T e fj = 3.5 MeV. The main effect was to increase the statistics of CC and ES and, above all, of 
NC events, and to increase significantly the precision on both the total 8 B neutrino flux and the 
neutrino mixing parameters. The value for the total 8 B neutrino flux extracted from neutral current 
was $atc = 5.14+q' 2 q x 10 6 cm _2 s _1 , where the error, obtained by summing in quadrature the statistic 
and systematic contributions, was reduced by more than a factor of two with respect to previous 
publications. For SNO data alone (LETA plus SNO III) the best fit point moved to the LOW region 
of parameter space, but the significance level was very similar to the one of the usual LMA solution. A 
global fit, including all the solar and the KamLAND data, essentially confirmed, instead, the previous 
results [13] for Am 2 2 and it made possible a further improvement in the angle determination, giving, 
in a 2 flavor analysis, tan 2 #i 2 = 0.457lo;o28- 

In the last five years also the Super-Kamiokande collaboration presented new analyses, including 
the data of the different working phases of this experiment: Super-Kamiokande II [68 1 (from December 
2002 to October 2005) and Super-Kamiokande III (from July 2006 to August 2008) [69] • Due to 
the 2001 accident, which damaged some of the photomultiplier tubes, the detector sensitivity was 
reduced with respect to SK-I and therefore it was important to improve the methods adopted for data 
collection (particularly for vertex event reconstruction, angular resolution and background reduction) 
and analysis. In this way, during the 548 days of SK-III a 2.1% systematic uncertainty on the total 
flux (corresponding roughly to two thirds of the SK-I value) was reached. The second and third Super- 
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Kamiokande phases essentially confirmed the SK-I results, for what concerns the absence of significant 
spectral distortion, the total 8 B measured flux and the day-night asymmetry. 

Since September 2008, Super-Kamiokande is running with modernized data acquisition system 
(DAQ) and electronics, which allow a wider dynamic range in the measured charge and is read out via 
Ethernet. This phase of the experiment is denoted as Super-Kamiokande- IV [70J. Thanks to the fast 
DAQ every hit can be recorded and the resulting data stream analyzed by an online computer system 
that finds timing coincidences which are saved as triggers. As a consequence, Super-Kamiokande's en- 
ergy threshold is now only limited by computing speed and the event reconstruction. The present event 
reconstruction is able to reconstruct electrons with a total energy of 3 MeV or more. The computing 
speed limits the energy threshold to 4.2 MeV which is just below the threshold of Super-Kamiokande-I 
and III (4.5 MeV). The same water flow techniques developed during Super-Kamiokande-III result in 
an observed solar neutrino elastic scattering peak between 4 and 4.5 MeV total recoil electron energy. 
Special techniques are developed to discriminate the signal from the background, taking advantage 
from the fact that the background is mainly due to (3 emission from 214 Bi and it is characterized 
by a larger Coulomb multiple scattering. This makes possible a reduction of about 10 — 15% of the 
statistical uncertainty and this method can also be applied to previous phases of the experiment. The 
additional systematic uncertainty of this method is under investigation. 

14.5.2 The impact of KamLAND results on solar neutrino physics 



Even if it is based on the analysis of a reactor antineutrino beam, the KamLAND experiment 
played a fundamental role in the solution of the long standing solar neutrino puzzle. In fact, the 
first KamLAND data [16] were determinant, in conjunction with the previous solar neutrino experi- 
ments (and mainly with SNO) and assuming CPT invariance, to prove the validity of the oscillation 
hypothesis and to select the LMA solution as the correct one. 

Between March 2002 and January 2004 a new set of data were collected and the KamLAND col- 
laboration performed a study including also a re-analysis of the previous data. During the 2002-2004 
campaign important upgrades were done both on the central detector (increasing the photocatode 
coverage and improving the energy resolution) and in the analysis techniques (reduction of the back- 
ground with better techniques in the event selection cuts based on the time, position and geometry 
of the events). The number of antineutrino events above 2.6 MeV expected in absence of antineutrino 
disappearance was 365.2 ± 23.7(syst) and the 258 observed events corresponded to a i7 e survival prob- 
ability equal to 0.658 ± 0.044(stat) ± 0.047(syst). The energy spectrum analysis was in disagreement 
with the no oscilation hypothesis at 99.6% statistical significance. In [19] the KamLAND collabora- 
tion, looking at the Lq/E spectrum dependence (where Lq is the source-detector distance and E the 
v e energy), performed also an interesting study of other alternative hypotheses (like decoherence and 
decay) for neutrino disappearance. The oscillation hypothesis offered by far the best explanation of 
the spectrum shape, as one can see from Fig. (14.3). 

As shown in Fig. ( 14.4 A.), the best fit obtained from the data analysis was in the so-called LMAI 
region (with values of Amj 2 around 8 • 10 -5 eV 2 ) and the alternative solution at higher Am 2 2 (around 
2 • 10 _4 eV 2 ) was strongly disfavoured, at 98% C.L., mainly due to the spectrum distorsions. The 
KamLAND data alone were not sufficient to solve completely the ambiguity on the mixing angle values 
and to exclude maximal mixing. However, including in the analysis also the results coming from solar 
neutrino experiments, the allowed values of the angle were significantly restricted (see Fig.( 14.4^ 3)) 
and the two flavor combined analysis gave Amf 2 = 7.91q5 • 10 -5 eV 2 , tan 2 #i2 = 0.40^q'qy at a 1 a 
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Figure 14.3: Ratio of the observed v e spectrum to the expectation for no-oscillation versus Lo/E. The 
curves show the expectation for the best-fit oscillation, best-fit decay and best-fit decoherence models, 
taking into account the individual time-dependent flux variations of all reactors and detector effects. 
Taken from [19J. 
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Figure 14.4: (A) Allowed region of the neutrino oscillation parameter from KamLAND anti-neutrino 
data (colored regions) and solar neutrino experiments (lines) [12] . (B) Result of a combined two- 
neutrino oscillation analysis of KamLAND and the observed solar neutrino fluxes under the assumption 
of CPT invariance. Taken from [19J. 

The next KamLAND analysis [20] included also, in addition to the one of [161 119] - the new data 
collected up to May 2007. The increase in data collection was significant (also thanks to the enlarging 
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the radius of the fiducial volume from 5.5 to 6 m) and there was a reduction of systematic uncertainties, 
in the number of target protons and the background. The total uncertainty on Am|i was around 2%, 
mainly due to the distortion of the energy scale in the detector. The total uncertainty, 4.1%, on 
the expected event rate was due to different sources (above all the definition of the detector fiducial 
volume and energy threshold, the v e spectra and the reactor power) and it affected primarily the 
mixing angle determination. The different background sources were studied and reduced further. The 
most important one was the 13 C(ce, n) 16 reaction, made possible by the a decay of 210 Po (a daughter 
of 222 Rn) introduced in the liquid scintillator during the construction, which produces neutrons with 
energies up to 7.3 MeV. 

The results of the statistical analysis are reported in Fig. ( 14.5) ), taken from [20]. The allowed 
oscillation parameter values were Am 2 x = 7.58 +ai4 ' " ■ 

,2/ 



_u 13 (stat) 

,2 



+0.15 



o.i5 ( s y st ) 



10 5 eV 2 for the mass eigenval- 
ues and tan 2 #i2 = 0.56^QQ7(stat)+QQg(syst), for tan 2 #i2 < 1 and the no oscillation hypothesis was 
excluded at 5 a. The extension to the three neutrino oscillation analysis had the main effect to enlarge 
the uncertainty on #12, leaving Am 2 2 substantially unchanged. Figure (14.5), taken from [20], shows 
that the effect of the inclusion in the analysis of the data from SNO [62J and previous solar neutrino 
experiments was essentially to reduce the interval of allowed 612 values and also to move the best fit 
point towards slightly lower values of the mixing angle. 
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Figure 14.5: Allowed region for neutrino oscillation parameters from KamLAND and solar neutrino 
experiments. The side-panels show the Ax 2 -profiles for KamLAND (dashed) and solar experiments 
(dotted) individually, as well as the combination of the two (solid). Taken from [20] . 



Figure (14.6) (taken from [20]) illustrates, instead, the v e survival probability, as a function of 



the ratio Lq/E between the average baseline and the antineutrino energies. One can notice that the 
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observed spectrum (after subtraction of background and geo-neutrino signals), reproduces correctly 
the general shape of the expected oscillation cycle, with a slight excess of low energy antineutrinos, 
that could be interpeted as geo-neutrinos. 



1 

| 0.8 
o 

£ 0.6 

> 

| 0.4 
0.2 



" 20 30 40 50 60 70 80 90 100 
L,/E v (km/MeV) 

Figure 14.6: Ratio of the background and geo-neutrino-subtracted 9 e spectrum to the expectation 
for no-oscillation as a function of Lq/E. Lq is the effective baseline taken as a flux- weighted average 
(Lq = 180 km). The energy bins are equal probability bins of the best-fit including all backgrounds. 
The histogram and curve show the expectation accounting for the distances to the individual reactors, 
time-dependent flux variations and efficiencies. The error bars are statistical only and do not include, 
for example, correlated systematic uncertainties in the energy scale. Taken from |20j . 




14.5.3 Toward the sub-MeV analysis: the Borexino detector and its measurements 

In the last decade significant steps forward have been done in the knowledge of solar neutrino 
properties, thanks mainly to the results obtained by the kiloton scale Cerenkov detectors (SK and 
SNO) and by advent of the reactor neutrino experiment KamLAND. However, these experiments 
investigated only the energy part of solar neutrino spectrum above 5 MeV, which represents a small 
fraction of the full spectrum. The single components of the neutrino spectrum cannot be determined 
by such techniques at low energies and, therefore, up to the last four years, low energy neutrinos had 
been observed only via radiochemical methods. A significant change took place with the advent of 
Borexino, a real time experiment which opened the way to the investigation of the sub-MeV region 
and isolated for the first time the neutrinos corresponding to the monochromatic berillium line. 

The Borexino detector 

Borexino is an ultra-high radiopure large volume liquid scintillator detector (using pseudocumene 
-PC-Q as aromatic scintillation solvent, and PPCj^] as solute at a concentration of 1.5 g/1) located 

7 1 ,2,4-trimethylbenzene 
8 2,5-diphenyloxazole 
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underground at the italian Gran Sasso National Laboratories (LNGS), under about 1400 m of rock 
(3800 mwe) |71j. The employment of a liquid scintillator as target mass assures a light production 
sufficient to observe low energy neutrino events via elastic scattering by electrons. This reaction is 
sensitive to all neutrino flavors, through the neutral current interaction, but the cross section for v e 
is larger than v„ and v T by a factor of 5-6, due the combination of charged and neutral currents. The 
main goal of Borexino is the measurement of the mono-energetic (0.862 MeV) 7 Be neutrinos, which 
have the basic signature of the Compton-like edge of the recoil electrons at 665 keV (see Fig. 14.7). 



■ All solar neutrinos 
- 7 Be neutrinos 




0.8 1 1.2 1.4 
Energy [MeV] 



Figure 14.7: Neutrino spectra expected in Borexino (accounting for the detector's energy resolution). 
The upper line represents the neutrino signal rate in Borexino according to the most recent predictions 
of the Standard Solar Model [72] including neutrino oscillations with the LMA-MSW parameters. 
The lower line illustrates the contribution due to Be neutrinos. The pp neutrinos contribute to the 
spectrum below 0.3 MeV and the edge at 1.2 MeV is due to pep neutrinos (from |73j). 



The high light yield typical of a liquid scintillator makes it possible to reach a low energy threshold, 
a good energy resolution of about 5% at 1 MeV and a pulse shape discrimination between a and f3 
decays. On the other hand, no directionality is possible and it is also not possible to distinguish 
neutrino scattered electrons from electrons due to natural radioactivity. For this reason, an extremely 
low level of radioactive contamination is compulsory and this has been one of the main tasks and 
technological achievements of the experiment. The background due to the presence of j3 decay of 14 C 
{fi en d-point 156 keV), intrinsic to the scintillator, limits neutrino observation to energies above 200 keV. 
Techniques for the scintillator purification are based mainly on methods developed and tested in earlier 
studies with the Counting Test Facility (CTF), a 4-ton prototype of Borexino which demonstrated for 
the first time the feasibility of achieving the low backgrounds needed to detect solar neutrinos in a 
large scale scintillator [741 175} 176]. For Borexino, a larger purification plant was developed similar to 
the CTF system, but with several improved features including the use of high vacuum and precision 
cleaning techniques. 

The design of Borexino is based on the principle of graded shielding (onion-like structure - see 
Fig. [148]). 

The scintillator (~ 300 tons) is contained in a thin nylon Inner Vessel (IV), of radius 4.25 m, 
at the center of a set of concentric shells of increasing radiopurity and it is surrounded by an outer 
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Figure 14.8: Schematic view of the Borexino detector. 



vessel (OV), filled with PC and 5.0 g/1 DMP |^J a material which is able to quench the residual 
scintillation of PC and acts as a passive shield against radon and other background contaminations 
originating from the external parts. A third more external vessel is composed of a stainless steel sphere 
(SSS), enclosing the passive shield (PC-DMP), and the entire detector is contained in a dome-shape 
structure 16.9 m high with a radius of 9 m, filled with ultra-pure water, denominated Water Tank 
(WT). The scintillation light is recorded by 2212 8-inches photomultipliers distributed on the inner 
part of the SSS [77J EE]; 1828 of them are equipped with aluminum light concentrators designed to 
increase the light collection efficiency [79] ■ Cerenkov light and residual background scintillation in 
the buffer are thus reduced. The others 384 photomultipliers without concentrators are used to study 
this background and to identify muons that cross the buffer and not the Inner Vessel. The Water 
Tank is equipped with 208 8-inches photomultipliers and acts as a Cerenkov muon detector. Although 
the muon flux is reduced by six order of magnitude by the 3800 m.w.e. depth of the Gran Sasso 
Laboratory, is still significant (1.1 fi mT 2 hr 1 ). An additional reduction, of the order of about 10 4 , 
has been necessary; for more details see Ref. [80]. 

In order to remove contaminants from dust (U, Th, K), air ( 39 Ar, 85 Kr) and cosmogenically 
produced isotopes ( 7 Be), different purification techniques were applied, such as distillation, water 
extraction, nitrogen stripping and ultra-fine filtration. The pseudocumene was distilled in-line during 
the detector filling at 80 mbar and at a temperature of about 90-95 °C. Distilled pseudocumene was 
stripped in a 8 m-high (15 cm in diameter) packed column with specially prepared ultra-low Ar/Kr 
nitrogen (0.005 ppm Ar and 0.06 ppt Kr, see Ref. |81j). Position reconstruction of the events, as 
obtained from the photomultipliers timing data via a time-of-flight algorithm, allowed to define a 
fiducial spherical volume, corresponding approximately to 1/3 (i.e. about 100 tons) of the scintillator 
volume in order to reject external 7 background. The others 2/3 of the scintillator act as an active 
shield. 

9 dimcthylphthalate 
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The measurement of the Be line 



The Borexino collaboration started taking data in May 2007 and after only 3 months (47.4 live 
days) it was able to extract the 7 Be signal from the background. The best value estimate for the rate 
was 47 ± 7(stat) ± 12 (syst) counts/(day • 100 ton), where the systematic error is mainly due to the 
fiducial mass determination [73J. An update of the 7 Be signal was reported after 9 months from an 
analysis of 192 live days (from May 16 th 2007 to April 12 th 2008), corresponding to 41.3 ton-yr fiducial 
exposure to solar neutrinos. 

The severe cuts that had to be passed by the events in order to be selected and enter the anal- 
ysis were mainly designed to avoid pile up of multiple events, reject the events originated by muons 
and their daughters and the ones due to radon daughters preceding the a — (3 Bi-Po delayed coinci- 
dences. Moreover, severe cuts (radial and based on the z-coordinates) were finalized to reduce the 
external 7 background. The remaining fiducial mass was of 78.5 tons. Important background sources 
were the fast coincidence decays from the 238 U chain (contamination level of (1.6 ± 0.1) 10 -17 g/g) 
and the 232 Th chain (contamination level of (6.8 ± 1.5) 10~ 18 g/g) and the 85 Kr contained in the 
scintillator that produces the rare decay sequence 85 Kr — > 85m Rb + e + + v e , 85m Rb — > 85 Rb + 7. 
The total estimated systematic error was 8.5% [65j . mainly determined by two sources, introduc- 
ing an uncertainty of 6% each: the total uncertainty on the fiducial mass and the one on the re- 
sponse function. The best value for the interaction rate of the 0.862 MeV 7 Be solar neutrinos was 
49 ± 3(stat) ± 4(syst) counts/ (day- 100 ton). This result excludes at the 4 a C.L. the no oscillation 
hypothesis for 7 Be solar neutrinos, which in the high metallicity SSM |82| [72] would imply 74 ± 4 
counts/ (day- 100 ton). The Borexino result is, instead, in very good agreement with the predictions of 
the LMA oscillation solution: 48 ± 4 counts/ (day- 100 ton). 



In order to reduce the systematic uncertainties and to tune the reconstruction algorithm and Monte 
Carlo simulations, a calibration campaign was performed in 2009 introducing inside the Borexino 
detector several internal radiosources a's, /3's, 7's, and neutrons, at different energies and in hundreds of 
different positions, which were determined with a precision better than 2 cm. The previous systematic 
error on 7 Be solar neutrino flux was estimated to be |65j at the level of 6% for both the fiducial 
volume and the energy scale. In the calibration campaign, the detector energy response was studied 



with eight 7 sources and Am-Be neutron source and comparing the calibration data and Monte 



Carlo simulations at different energies within the solar neutrinos energy region. The energy scale 
uncertainty, obtained with these studies, was determined to be less than 1.5%. 

The inaccuracy of the position (reduced by means of studies with a and /3 events) was less than 3 
cm, equivalent to a systematic error of 1.3% for the overall fiducial volume in the 7 Be solar neutrino 
energy region. The analyzed data set run from May 2007 to May 2010, with a fiducial exposure 
equivalent to 153.6 ton-year. In order to extract the 7 Be solar neutrino signal, the spectral fit was 
applied assuming all the intrinsic background components such as 85 Kr, 210 Bi, 14 C, n C. The Be solar 
neutrino rate was evaluated to be 46.0 ± 1.5(stat) ± 1.3(syst) counts/day- 100 ton [83]. Thanks to the 
calibration campaign, the systematic error was reduced to 2.7% and the total uncertainty to 4.3%. 



14.5.4 The pep and CNO neutrinos measurement in Borexino 



"When thermal neutrons are captured by protons a 2.2 MeV 7-ray is generated. 
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Figure 14.9: Examples of fitted spectra; the fit results in the legends have units [counts/(day-100ton)]. 
Left panel: A Monte Carlo based fit over the energy region 270-1600 keV to a spectrum from which 
some, but not all, of the a events have been removed using a PSA cut, and in which the event energies 
were estimated using the number of photons detected by the PMT array. Right panel: An analytic 
fit over the 290-1270 keV energy region to a spectrum obtained with statistical a subtraction and in 
which the event energies were estimated using the total charge collected by the PMT array. In all 
cases the fitted event rates refer to the total rate of each species, independently from the fit energy 
window (from |83j). 



In the SSM, due to the solar luminosity constraint and their intimate link to the pp neutrinos 
3"B"] . the mono-energetic 1.44 MeV pep neutrinos have one of the smallest uncertainties (1.2%) 
For this reason, after the pp neutrinos, they constitute the ideal probe to test SSM hypotheses. 
On the other hand, the detection of neutrinos within the CNO-bicycle is central to probe the solar 
core metallicity and contribute in this way to the solution of the solar metallicity problem |84t 148] . 
Also, they are believed to fuel massive stars with mass greater than ~ 1.2 M© during main sequence 
evolution and also stars with lower masses in more advanced stages of evolution. The energy spectrum 
of neutrinos from the CNO-bicycle is the result of three continuous spectra with end point energies of 
1.19 MeV ( 13 N), 1.73 MeV ( 15 0) and 1.74 MeV ( 17 F). Despite their relevance, until 2011, no pep and 
CNO neutrinos had been detected directly. 

The electron recoil energy spectrum from pep neutrino interactions in Borexino is a Compton-like 
shoulder with end point of 1.22 MeV, as one can see from Fig. (14.10), showing the pep and CNO 
contribution in Borexino . 

As already mentioned, very low background levels \73\ [65] are required to detect 7 Be neutrinos; 
the detection of pep and CNO neutrinos is even more challenging, as their expected interaction rates 
are ~10 times lower. The expected rate is on the order of a few counts per day in a 100 ton target. 
To detect pep and CNO neutrinos the Borexino Collaboration adopted a novel analysis procedure to 
suppress the dominant background in the 1-2 MeV energy range, due to the cosmogenic /3 + -emitter 
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Energy spectrum simulation 
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Figure 14.10: The neutrino-induced electron recoil spectra expected in Borexino. The total rates are 
those predicted by the latest high-Z solar model [48J . The pep and CNO neutrinos recoil spectra with 
end points in the region 1.2-1.5 MeV are shown. Also the Be neutrinos (measured in [83J), with a 
count rate about 10 times larger, are shown for comparison. Note that the variable on the x axis is 
not directly the energy value. Taken from [89j. 



11 C produced within the scintillator by muon interactions with 12 C nuclei. The muon flux crossing 
the Borexino detector, ~4300/i/day, yields a n C production rate of ~27 counts/ (day- 100 ton). This 
background can be reduced by performing a space and time veto following coincidences between signals 
from the muons and the cosmogenic neutrons [85, 86j, discarding exposure that is more likely to contain 
n C due to the correlation between the parent muon, the neutrorp*"| and the subsequent n C decay (the 
Three- Fold Coincidence, TFC). The TFC technique is based on the reconstructed track of the muon 
and the reconstructed position of the neutron-capture 7-ray [88]. The criteria of rejection were applied 
to obtain the best compromise between n C rejection and preservation of fiducial exposure, resulting 
in a n C rate of (2.5±0.3) count per day, (9±1)% of the original rate, while preserving 48.5% of the 
initial exposure. 



Figure (14.11) shows the resulting spectrum obtained with data collected between January 2008 



and May 2010, corresponding to a fiducial exposure of 20409 ton-day [58]. Despite the TFC veto, the 
number of 11 C surviving events still constituted a significant background. 

To discriminate 11 C /3 + decays from neutrino-induced e~ recoils and (3~ decays the pulse shape 
differences between e~ and e + interactions in organic liquid scintillators |90t l§T] were exploited. In 
fact a small difference in the time distribution of the scintillation signal arises from the finite lifetime 
of ortho-positronium as well as from the presence of annihilation 7-rays, which present a distributed, 
multi-site event topology and a larger average ionization density than electron interactions. The 



11 In 95% of the cases at least one free neutron is spalled in the 11 C production process 
scintillator with a mean time of 255 fj,s [88] . 



7\, and then captured in the 
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Figure 14.11: Energy spectra of the events in the FV before and after the TFC veto is applied. The 
solid and dashed blue lines show the data and estimated n C rate before any veto is applied. The solid 
black line shows the data after the procedure, in which the n C contribution (dashed) has been greatly 
suppressed. The next largest background, 210 Bi, and the electron recoil spectra of the best estimate 
of the pep neutrino rate and of the upper limit of CNO neutrino rate are shown for reference. Rate 
values in the legend are quoted in counts/(day • 100 ton) from [58] . 



Borexino Collaboration employed an optimized pulse shape parameter using a boosted-decision-tree 
algorithm [92], trained with a TFC-selected set of n C events (e + ) and 214 Bi events (e~) selected by 
the fast 214 Bi- 214 Po a-/3 decay sequence. In a work published in 2012 |58j the Borexino Collabora- 
tion presented the results of an analysis based on a binned likelihood multivariate fit performed on 
the energy, pulse shape, and spatial distributions of selected scintillation events whose reconstructed 
position is within the fiducial volum^J 

The energy spectra and spatial distribution of the external 7-ray backgrounds have been obtained 
from a full, Geant4-based Monte Carlo simulation, and validated with calibration data from a high- 
activity 228 Th source [93J deployed in the outermost buffer region, outside the active volume, a events 
were removed from the energy spectrum by the statistical subtraction method [73]. In the energy region 
of interest of the fit procedure all background species whose rates were estimated to be less than 5% of 
the predicted rate from pep neutrinos have been excluded. All rates were constrained to positive values 



and thirteen species were left free in the fit 13 The rate of the radon daughter 214 Pb was fixed using 
the measured rate of 214 Bi- 214 Po delayed coincidence events. The contribution from pp solar neutrinos 
was fixed to the SSM assuming MSW-LMA with tan 2 12 =0A7±%%1, Am 2 2 =(7.6±0.2)-10- 5 eV 2 
and the contribution from 8 B neutrinos to the rate from the measured flux [661 195] . 



12 less than 2.8 m from the detector center and with a vertical position relative to the detector center between -1.8 m 
and 2.2 m. 

13 electron recoils from 7 Be, pep, and CNO solar neutrinos, internal radioactive backgrounds 210 Bi, 11 C, 10 C, 6 He, 40 K, 
85 Kr, and 234m Pa, and external 7-rays from 208 T1, 214 Bi, and 40 K. 
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In Table 14.6 the results for the pep and CNO neutrino interaction rates are shown. The ab- 
sence of a pep neutrino signal was rejected at 98% C.L. Concerning the CNO neutrinos flux, its 
electron-recoil spectrum is similar to the spectral shape of 210 Bi, but the last one is about 10 times 
greater; therefore it has only been possible to provide an upper limit on the CNO neutrino inter- 
action rate. The 95% C.L. limit reported in Table 14.6 has been obtained from a likelihood ratio 
test with the pp neutrino rate fixed to the SSM prediction [38] under the assumption of MSW-LMA, 
(2.80±0.04) counts/ (day- 100 ton). 



V 


Interaction rate 


Solar- v flux 


Data/SSM 




[counts/ (day-100ton)] 


[lO^m^s- 1 ] 


ratio 


pep 


3.1 ± 0.6 sta t± 0.3 syst 


1.6 ±0.3 


1.1 ±0.2 


CNO 


< 7.9 (< 7.1 sta tonly) 


< 7.7 


< 1.5 



Table 14.6: Best estimates for the pep and CNO solar neutrino interaction rates. For the results in 
the last two columns both statistical and systematic uncertainties are considered. Total fluxes have 
been obtained assuming MSW-LMA and using the scattering cross-sections from |96[ l94"t [97] and a 
scintillator e~ density of (3.307±0.003)-10 29 ton" 1 . The last column gives the ratio between our 
measurement and the high-Z (GS98) SSM [38]. Table taken from [58] . 



14.6 Phenomenological analysis 



14.6.1 Status of the determination of the mixing parameters in a 3 flavor analysis 



Recently, the SNO collaboration performed a combined analysis of all the three working phases of the 
experiment [95] based on a fit to Monte Carlo derived probability density functions (PDFs) for each of 
the possible signals and backgrounds, and also introduced a new way to parametrize the 8 B neutrino 
signal. Figure (14.12), reporting the results of the two flavour (with the assumption #13 = 0) SNO 
only analysis, shows the further improvement in the mixing parameters accuracy, but, at the same 
time, it confirms that the SNO results alone would not be sufficient to completely exclude the LOW 
solution. 

This ambiguity was definitely removed, as shown in Figure (14.13), by including in the analysis 
the results of all previous solar neutrino experiments |98|. 199} 164"! 168] [69] , the 7 Be solar neutrino rate 
measured by Borexino [83], the 8 B neutrino spectra [100J and the KamLAND datsP] |21j . 



The higher values of Am 2 2 in the LMA region were excluded, together with the full LOW solution, 
thanks mainly to the large discrimination power of KamLAND. This experiment, however, did not 
contribute significantly to improve the mixing angle determination and the accuracy on this parameter 
remained quite high. The results of the two flavor analysis are reported in Table (14.7) (taken from 

m\). 



14 The KamLAND data were obtained in a completely independent experiment and, therefore, the corresponding \ 2 
values, as functions of the mixing parameters, were directly summed to the x 2 values computed by direct solar neutrino 
analysis. 
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Figure 14.12: Two-flavor neutrino oscillation analysis contour using only SNO data (taken from [95]). 




Figure 14.13: Two-flavor neutrino oscillation analysis contour using both solar neutrino and Kam- 
LAND results (taken from |95j). 



The slight tension between the solar neutrino experiments and KamLAND was significantly reduced 
by extending the analysis to the 3 flavor oscillation case as shown in Figure (14.14), from which it is 
clear that the best global fit is obtained for values of #13 different from zero. 

A detailed analysis of the \ 2 behavior proved also that the combination of solar experiments and 
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Analysis tan 2 8 12 Am^eV 2 ] X 2 / NDF 

SNO only (LMA) 5.62+}-^ x 10~ 5 1.39/3 

SNO only (LOW) 0.427±g;g|| l-35±°;^ x 10~ 7 1.41/3 

Solar 0.427±g;g|| 5.13±J;|| x 10~ 5 108.07/129 

Solar +KamLAND 0.427±g;g|X 7.46±g;gg x 10~ 5 

Table 14.7: Best-fit neutrino oscillation parameters from a two-flavor neutrino oscillation analysis. 
Uncertainties listed are 1 a after the x 2 was minimized with respect to all other parameters (taken 
from |95|). 




Figure 14.14: Three- flavor neutrino oscillation analysis contour using both solar neutrino and Kam- 
LAND results. Taken from [95j . 



KamLAND enables to improve significantly the discriminating power on the #13 mixing parameter 
(see Figure |14.15| and Table |14.8[ ) . 

The indication in favor of $13 being different from zero was in agreement with the recent results 
from the long-baseline experiments T2K |101j and MINOS |102| . and with the combined analysis 
performed in |103j . including also the atmospheric neutrino and the CHOOZ |104j data. Moreover 
the validity of this hint has been corroborated by the data obtained this year by the short baseline 
neutrino reactor experiments |105] 11061 1107] . which established that 6*13 > at about 5a (and even 
more in the Daya Bay case [108J). These experiments found values of sin 2 #13 centered between 0.020 
and 0.030; very promising results for future experiments looking for leptonic CP violation |109j . The 
impact and the possible consequences of these recent results have been discussed, among the others, in 
the following papers f!09l IllOl 1111] . The different accuracy that can be reached in the determination 
of the mixing angle between the first and third generation, according to the different kind of neutrino 
experiments included in the analysis, is represented in Figure (14.16). 

The combined analysis of the different SNO phases was also very useful to obtain a precise de- 
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Figure 14.15: Projections of the three-flavor neutrino oscillation parameters. The horizontal lines 
represent the A% 2 for a particular confidence level. Taken from [95] • 
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Table 14.8: Best-fit neutrino oscillation parameters from a three-flavor neutrino oscillation analysis. 
Uncertainties listed are ±1 a after the x 2 was minimized with respect to all other parameters. The 
global analysis includes Solar+KL+ATM+LBL+CHOOZ. 

termination of the 8 B solar neutrino flux, <&8 B = 5.25 ± 0.16(stat)^Q'J3(syst) x 10 6 cm~ 2 s~ 1 , with an 
important reduction of the systematic uncertainty. This result was consistent with, but more precise 
than, both the high-Z BPS09(GS), $ = (5.88 ± 0.65) x lO^m^s" 1 , and low-Z BPS09(AGSS09), 
$ = (4.85 ± 0.58) x 10 6 cm _2 s _1 , solar model predictions |112| . 

The combination of the LETA analysis by the SNO collaboration [66J and of the Borexino measure- 
ments |100| made possible a detailed study of the low energy part of the 8 B solar neutrino spectrum. 
Even if characterized by a larger uncertainty (mainly due to a more limited statistics), Borexino data 
confirm the LETA indication of low energy data points lower than the theoretical expections based 



on matter enhanced oscillation and solar models as shown in Figure 14.17| (taken from [100J). These 
results agreed also with the Super-Kamiokande observation [64J of flat spectrum, consistent with the 
undistorted spectrum hypothesis. The emergence of this slight tension between theory and experi- 
ments seems to indicate the presence of new subdominant effects and also suggests the possibility of 
non-standard neutrino interactions (like those studied in |113j ) or the mixing with a very light sterile 
neutrino |114j . Future solar neutrino experiments, like SNO+, could shed more light on this subject, 
by performing precision measurements of lower energies solar neutrinos (like the pep neutrinos) . 
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Figure 14.16: Projection over sin 2 #13 combining the projections obtained by analyzing data from 
all neutrino sources. The data from atmospheric, short-baseline experiments and long-baseline ex- 
periments (ATM+LBL+CHOOZ) was determined from Figure 2 (left panel) in [103] which already 
includes the latest T2K [101J and MINOS [TD2] results. 




14.6.2 Free flux analyses 

The increasing data of solar neutrinos allow to independently test the astrophysics of the solar interior 
and the physics of neutrino propagation. The analysis discussed in previous sections can be modified 
by also varying the solar neutrino fluxes in order to accommodate all neutrino data, while all the 
functional dependences are maintained as predicted by the standard model dependences. A key step 
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in this kind of analysis is the imposition of the luminosity constraint |115| 1116] . which implements 
in a global way for the Sun the constraint of conservation of energy for nuclear fusion among light 
elements. Each neutrino flux is associated with a specific amount of energy released to the star and 
therefore a particular linear combination of the solar neutrino fluxes is equal to the solar luminosity 
(in appropriate units). One can write the luminosity constraint as 

where L is the solar luminosity measured at the earth's surface, 1 A.U. is the average earth-sun 
distance, and the coefficient oti is the amount of energy provided to the star by nuclear fusion reactions 
associated with each of the important solar neutrino fluxes, <3?j. The coefficients 014 are calculated 
accurately in ref. [116J. 

The model independent determination of the solar neutrino fluxes |118| 1117] shows that present 
solar neutrino data leads to accurate results for four fluxes and also the correlations between them. 
This information allows for a consistent global comparison of SSM fluxes with the inferred fluxes by 
neutrino data. Present data leads to the values for the inferred solar neutrino fluxes reported in the 



fourth column (labelled as "Solar") of Table 14.4 in Section 14.3 The precision of the 7 Be and 8 B 



neutrino fluxes is driven by the Borexino and SNO (SK) neutrino experiments, while the precision 
of the pp and pep neutrino fluxes mainly comes by the imposition of the luminosity constraint. The 
neutrino data directly demonstrates that the Sun shines by the pp chain. The CNO cycle only 
contributes to the total luminosity at the percent level. 

The reader may wonder how much these inferences are affected by the luminosity constraint. The 
idea that the Sun shines because of nuclear fusion reactions can be tested accurately by comparing 
the observed photon luminosity of the Sun with the luminosity inferred from measurements of solar 
neutrino fluxes. Moreover, this same comparison will test a basic result of the standard solar model, 
namely, that the Sun is in a quasi-steady state in which the current energy generation in the interior 
equals the current luminosity at the solar surface. The free flux analysis, without imposing luminosity 
constraint, permits an estimation of the solar luminosity inferred by neutrino data, which agrees with 
the directly measured one within 15 % (1 a). 

14.7 Future solar neutrino experiments 

14.7.1 The near future: improvement of pep measurements and CNO detection 

In the last decades the intensive study of 8 B and, more recently, 7 Be solar neutrinos made possible 
fundamental steps forward in the solution of the solar neutrino puzzle and the determination of 
the neutrino mixing parameters. Nevertheless, many key features of the oscillation models (like the 
transition between the vacuum dominated sub-MeV region and the spectral region between 1 and 3 
MeV, where matter effects become relevant) still have to be tested or verified with better accuracy 



and precision (see Figure 14.18, taken from [119J). 



The apparent partial deficit of events in the low energy part of the B spectrum suggested the 



introduction of new theoretical models (as discussed in section 14.6). Also for these reasons, the 
experimental efforts in the last years focused on the detection of neutrinos of ever decreasing energies, 
to fully confirm the validity of the MSW-LMA solution and verify the fluxes predicted by SSMs, 
discriminating between different version of these models. The fluxes of the medium and high energy 
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Figure 14.18: The u e survival probability is represented as a function of neutrino energy. The gray 
band represent the MSW-LMA prediction. The higher survival probability region at low energies is 
where vacuum- dominated oscillations occur. As the neutrino energy increases, matter effects become 
important and the lower survival probability at high energies is due to matter-enhanced oscillations. 
The reported data correspond to solar nuetrino flux measurements performed by different experiments. 
Taken from f[TT9]). 



neutrinos of the pp chains ( 7 Be, 8 B and hep) are predicted with quite large uncertainties, mainly due 
to the uncertainties in nuclear cross sections and solar opacity (Table 14.5). The pp and pep fluxes, 
instead, are strongly correlated between themselves and their values are predicted with the highest 
precision because SSMs predict that pp chain reactions are responsible for more than 99% of the energy 
powering the Sun [59]. Therefore, the measurements of these components would be the most stringent 
test of the SSM. The tight correlation between pep and pp neutrinos is theoretically well established 
and, therefore, even in the pessimistic hypothesis that pp neutrinos could not be measured with the 
desired accuracy, a significant improvement in the pep neutrinos measurement with respect to data 
presently available would make possible to reduce significantly the 15% indetermination on the solar 
luminosity (see subsection 14.6.2) and to test indirectly the SSM's predictions that almost 100% of 
solar energy is produced by nuclear burning. 

As already mentioned, water Cerenkov detectors, which played a fundamental role in the solution 
of the solar neutrino problem, are characterized by a low photon yield |1201 1121| and therefore can 
detect only the higher part of the spectrum (hep and 8 B neutrinos with a threshold around 3.5 MeV). 
The radiochemical experiments [99, I122j are limited, instead, by their ability to measure only the 
integrated neutrino rate above the charge-current interaction threshold (down to 0.23 MeV for the 
Gallium experiments), without the possibility to discriminate between the different spectrum compo- 
nents. Therefore, an important contribution should come from the present and future organic liquid 
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scintillator detectors, planned to perform low energy solar neutrino spectroscopy. To reach this goal, 
they will take advantage from the high values of light yield (about 10 4 photons per MeV of deposited 
energy) and from the possibility to assemble very large masses of high purity material. The excellent 
levels of radiopurity, reached for instance at Borexino, and the typical geometry of these detectors 
(which are unsegmented and can be easily adapted to the definition of a fiducial volume) are funda- 
mental to reduce the impact of the background, that is so critical due to the feebleness of the low 
energy signal. 

In the near future significant contributions are expected from Borexino and SNO+ |123j exper- 
iments. Borexino has already proved its importance in this kind of analysis performing the first 
measurements of pep and CNO neutrinos (even if the level accuracy is not yet the desired one) and 
further reducing, with the purification campaign started since July 2010, the level of contamination 



from almost all of the main radioactive background sources 15 The purification efforts are still ongoing 
and should make possible a further improvement on the accuracy of the signal extraction. The SNO+ 
experiment, that should start taking data soon in the SNOLAB, should take advantage from the lo- 
cation (about two times deeper underground than the Gran Sasso laboratory), with the consequent 
lower muon flux and a strongly reduced n C rate. Moreover, thanks to the detector mass (about three 
times larger than in Borexino), it should be able to reach a higher counting rate. This could deter- 
mine a fundamental improvement at least in the case of the pep neutrino measurement, where a 5% 
uncertainty is expected, to a level that should make possible a significant test of the MSW transition 
region. 

In the more optimistic scenarios it may be also possible to attach the main problem of measuring 
lowest energy parts of the solar neutrino spectrum, that is the pp neutrinos and the 0.38 MeV Berillium 
line. In any case the presence in organic scintillators of an intrinsic 14 C background will make this 
very low energy measurements an extremely hard task and they may require the introduction of new 
techniques, like the ones we are going to describe in the next subsections. 

14.7.2 The far future: experimental challenges 

The challenge for all future experiments aimed at measuring the low energy part of solar neutrino 
spectrum is that of assembling experimental devices with low energy thresholds suitable to detect a 
low rate signal in a region characterized by different potential sources of radioactive background. This 
difficult experimental task is common also to the experiments looking for neutrinoless double /3 decay 
or for dark matter signals (search for signatures of WIMPs, a stable or long-lived weakly interacting 
elementary particle, produced in the early Universe, whose existence is predicted in extensions of 
the Standard Model). In fact, some of the solar neutrino experiments planned for the future are 
multipurpose experiments designed also for the other above-quoted topics. 

They are all characterized by a very large detector target mass and by the need to reach very high 
levels of radiopurity. The common feature is that of using scintillator detectors, but they differ for 
the chosen active scintillator material, which can vary from traditional organic scintillators (developed 
with the use of innovative technological devices) to new materials, like the noble gases. 

Noble liquid detectors: CLEAN and XMASS 



3 The main problem still surviving seems to be the reduction of 1 Pb. 
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One of the possible future frontiers is the idea to use scintillation detectors with liquid noble gases, 
like xenon, argon and neon. These materials have the advantage of being relatively inexpensive, easy to 
obtain and dense and it is not too difficult to build large homogeneous detectors of this kind; moreover, 
they can be quite easily purified, offer very high scintillation yields (about 30 — 40 photons/keV) and 
do not absorb their own scintillation light. 



i 
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Figure 14.19: Scheme of the CLEAN detector. Taken from |124| . 



A first example is offered by the CLEAN/DEAP family, a series of detectors based entirely on 
scintillation in liquid neon (LNe) and liquid argon (LAr). They have been realized using a scaleable 
technology in order to reach increasing sensitivities in the different prototypes realized and installed in 
the SNOLAB (Pico-CLEAN, Micro-CLEAN, DEAP-I, Mini-CLEAN and CLEAN/DEAP) with the 
aim to search for dark matter and to perform (through the analysis of elastic neutrino-electron and 
neutrino-nucleus scattering) a real time measurement of the pp solar neutrino flux. The final detector 
CLEAN (Cryogenic Low Energy Astrophysics with Noble gases) |124j (see figure 14.19) will be made 
by a stainless steel tank, of about 6 meters of diameter, filled with 100 tons of cryogenic liquid neon; 
only the central part of it, surrounded isotropically by a series of photomultipliers, will constitute the 
detector fiducial volume. An external tank of water, 10 metres wide and 12 metres high, will act as 
7-ray shielding, neutron shielding and muon veto. According to Monte Carlo simulations, there should 
be a production of 15000 photons/MeV and it should be possible to reach a 100% photon wavelength 
shifter efficiency and a statistical uncertainty on the pp measurements of the order of 1 %. 

A precise measurement of the pp component and of the ratio between pp and 7 Be fluxes would be 
essential to test the predictions of SSMs. A high accuracy on the pp neutrino flux would also make 
possible a better determination of the 812 mixing angle, which, complemented with the results from 
previous solar neutrino experiments and from KamLAND (essential for the Amf 2 measurement), 
would be fundamental to test the consistency of the LMA solution also in the region of transition 
between vacuum dominated and matter enhanced oscillations. Finally, CLEAN could in principle try 
to measure also the CNO neutrino flux, through the analysis of neutrino spectrum from 0.7 to 1.0 
MeV, with an estimated accuracy between 10 and 15%. 
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Figure 14.20: Schematic view of the full XMASS facilty (left) and a detail of the inner detector (right 
panel), from which one can see the particular configuration of the hexagonal photomultiplier tubes. 
Taken from [127J and [128] . 



An interesting alternative to the use of neon is offered by liquid xenon scintillator detectors |125| . 
which take advantage of the fact that among liquid rare gases xenon has the highest stopping power 
for penetrating radiation (thanks to its high atomic number, A ~ 131 and density, p = 3g/cm 3 ) and 
also the highest ionization and scintillation yield. The technological improvements of the last twenty 
to thirty years made possible significant improvements in the cooling and purification techniques of 
this kind of detectors and in the possibility of assembleing large mass detectors, of the order of some 
tons (like in the case of MEG [126] experiment, studying the \i — > erf decay). 



The XMASS experiment (see Figure 14.20) is a multipurpose low background and low energy 
threshold experiment that will use a large massive liquid xenon detector and has been designed to 
look for WIMPs (dark matter candidates), search for neutrinoless double (3 decay and study the pp 
and the 7 Be solar neutrinos. After two preliminary phases, during which smaller prototypes have 
been realized and installed in the Kamioka mine |127j . and the first data on double beta decay and 
dark matter have been taken, the full XMASS detector (that will measure also solar neutrinos) will 
have a total mass of 20 tons, with a fiducial volume of 10 tons. Special efforts are required mainly 
to lower the background, by reducing the radioactive contamination in the parts used for detector 
construction (with special attention to the photomultipliers and the copper material used for PMT 
holder), constructing a larger pure water active shield (for muons and mainly neutrons and 7 rays) 
and, above all, developing a distillation system for xenon in order to reduce the contamination by 
85 Kr, the major source of radioactive background inside the detector. 

Another interesting experimental project based on the noble gases liquid scintillator technique 
is that of DARWIN (DARk matter Wimp search with Noble liquids) [129J, which brings together 
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differen European and US research groups working on existing experiments and on the study for a 
future multi-ton scale LAr and LXe dark matter search facility in Europe. The main goal of the 
experiment is to look for a WIMP signal and to demonstrate its dark matter nature, taking advantage 
from the fact of performing the measurement with multiple different targets operating under similar 
conditions. In this way, it should be possible to estimate the dependence of the rate with the target 
material and, therefore, to better determine the WIMP candidate mass and to distinguish between 
spin independent and spin dependent couplings. The energy region of the nuclear recoil spectrum, 
below 200 keV, that should be investigated by this future experiment is of particular interest also for 
the study of the pp solar neutrinos and, in fact, the elastic scattering on electrons by the low energy 
component of the neutrino spectrum would be one of the main background sources for WIMP searches 



in liquid xenon detectors, as shown in Figure 14.21 




Energy [keV] 

Figure 14.21: Expected nuclear recoil spectrum from WIMP scatters in LXe for a spin-independent 
WIMP-nucleon cross section of 10 -47 cm 2 (red solid) and 10~ 48 cm 2 (red dashed) and a WIMP mass 
of 100GeV/c 2 , along with the differential energy spectrum for pp (blue) and 7 Be (cyan) neutrinos, 
and the electron recoil spectrum from the double beta decay of 136 Xe (green). Assumptions are: 
99.5% discrimination of electronic recoils, 50% acceptance of nuclear recoils, 80% flat analysis cuts 
acceptance. Taken from the second paper of |129j . 



DARWIN officially started in 2010; a technical design study should be ready in Spring 2013 and 
the start of the first physics run is expected by mid 2017. 

Multi kiloton scale liquid scintillators: example LENA 

The Borexino experiment demonstrated the great potential of the liquid-scintillator technique for 
the detection of low energy solar neutrinos. Thanks to this experience, a next-generation neutrino 
detector has been proposed: LENA (Low Energy Neutrino Astronomy) |130j . LENA is a multipurpose 
detector aiming to study supernova neutrinos, diffuse supernova neutrino background, proton decay, 
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atmospheric neutrinos, long-baseline neutrino beams, geoneutrinos and, last but not least, solar neu- 
trinos. The LENA project foresees a cylindrical detector with a diameter of 30 m and a length of 
about 100 m. Inside the detector is foreseen an internal part (with a diameter of about 26 m) con- 
taining about 50 kilotons of liquid scintillator, separated from a non-scintillating buffer region by a 
nylon barrier. Outside, a tank (made in steel or concrete) separates the inner detector from an outer 
water tank; it is used both for shielding and as an active muon veto. To collect the scintillation light, 
about 45,000 photomultipliers (with a diameter of 20 cm) are mounted to the internal walls of the 
detector. To increase the optically active area, the photomultipliers tubes are equipped with conic 
mirrors, the corresponding surface coverage is about 30%. Figure (14.22) shows a schematic overview 
of the current LENA design. 



height: 115 m, diameter: 50 m 
shielding from cosmic rays: -4,000 m.v 
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Figure 14.22: Schematical view of the LENA detector. From [130J 



Among the favored solvent for the liquid scintillator in LENA, the LAB (linear alkylbenzene) is 
currently the preferred one. It has a high light yield and large attenuation length and it has also the 
advantage of being a non-hazardous liquid. The attenuation lengths is on the order of 10 to 20 m (at 
a wavelength of 430 nm) and the photoelectron yield could be greater than 200 photoelectrons per 
MeV (with a scintillator mixture containing 2g/l PPO and 20 mg/1 bisMSB as wavelength shifters). 
Studies have been carried out to test the large-scale light transport and the differences in scintillator 
response for a, f3 and 7 particles. An alternative solvent option is the well studied PXE [131 j or a 
mixture of PXE and dodecane. 

As already pointed out, Borexino has splendidly demonstrated the potential of the detection tech- 
nique with liquid scintillator based detectors for solar neutrino detection. This technique offers the 
opportunity for a spectrally resolved measurement of the solar neutrino spectrum in the all energy 
range. 



Because the smaller ratio of surface to volume compared to the Borexino detected 16 ! in LENA it 



16 A smaller ratio of surface to volume decreases the chance that the scintillator is contaminated with radioimpuritics 
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Source 


EW [MeV] 


m M [kt] 


Rate [cpd] 


PP 


>0.25 


30 


40 


pep 


0.8-1.4 


30 


2.8xl0 2 


7 Be 


>0.25 


35 


l.OxlO 4 


8 B 


>2.8 


35 


79 


CNO 


0.8-1.4 


30 


1.9xl0 2 



Table 14.9: Expected solar neutrino rates in LENA (channel ue — > ev). The estimates are derived 
from the existing Borexino analyses [132J HOOj as well as expectation values for the respective energy 
windows (EW) of observation (133|. 1134] f!35j. The quoted fiducial masses, rn^, in LAB are based on 
a Monte Carlo simulation of the external 7-ray background in LENA. Table taken and adapted from 



is very likely to reach the excellent background conditions of Borexino. Monte Carlo simulations of 
the gamma background due to the uranium, thorium and potassium from the photomultipliers glass 
shows that a fiducial volume of the order of 30 ktons is achievable for solar neutrino studies; LENA 
will be able to address topics both in neutrino oscillations and in solar physics thanks to its unprece- 
dented statistics. A high statistics can be obtained in short times and in both Pyhsalmi and Frejus 
underground laboratories, where the detector could be hosted, where the cosmogenic background of 
11 C will be significantly lower than in Borexino. 

Monte Carlo simulations show that for pep, CNO and low-energy 8 B-z^s detection a fiducial mass 
of ~30 kton is necessary, while the fiducial mass for 7 Be-^s and high-energy (E > 5 MeV) 8 B-i/s could 
be enlarged to 35 kton or more. 



In Table 14.9 are reported the expected rates in 30 kton for the neutrinos emitted in the pp chain 
and the CNO-bicycle, using the most recent solar model predictions. This evaluation refers to a 
detection threshold set at about 250 keV. 

New techniques with organic scintillators: LENS 



The main goal of the Low Energy Neutrino Spectroscopy (LENS) detector is the real time mea- 
surement of solar neutrinos as a function of their energy, focusing, in particular, in the analysis of the 
lowest energy neutrinos coming from proton-proton fusion (i.e. the pp neutrinos), which represent the 
main contribution and the less known component of the pp-chain of fusion reactions inside the Sun. 

In order to make an energy spectrum measurement on low energy neutrinos, it is necessary to reach 
a low threshold for the charged current (CC) process and to be able to discriminate the background 
from radioactive decays. The CC process employed in LENS is the neutrino induced transition of 
115 In to an excited state of 115 Sn: 

v e + 115 In ^ 115 Sn* + e" (E = E„ - 114keV) . (14.14) 

115 Sn*(r = 4.76/is) ^ 115 Sn + 7 (498keV) + 7 (116keV) . (14.15) 

Thanks to that it is possible to detect low energy neutrinos with a threshold of 114 keV and 
measure their energy, following an idea that has been investigated since the 1970's [136J. 
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The primary interaction and the secondary cascade enable a triple coincidence, correlated in space 
and time. LENS employs as detection medium a liquid scintillator chemically doped with natural 
indium ( 115 In = 95.7%). In order to exploit the spatial correlation, the volume of the detector is 
segmented into cubic cells (7.5 cm) by clear foils (Teflon FEP) that have a lower index of refraction 
than the liquid scintillator. By internal reflection, the scintillation light produced in a cell is channeled 
in the directions of the 6 cell faces. The collected channeled light is read-out at the edge of the detector 
by photomultiplier tubes. 

LENS should be able to determine the low energy solar neutrino fluxes with an accuracy < 4%, 
testing neutrino and solar physics with a global precision better than the present one and also looking 
for any inconsistency in the LMA conversion mechanism [137J. 

14.8 Open questions in solar neutrino physics 



14.8.1 The metallicity problem 

The solar abundance or solar metallicity problem has been around for some time now. In analogy 
with the solar neutrino problem, there have been attempts (although in most cases, it is fair to say, of 
somewhat less radical nature) to solve it by introducing modifications to the input physics of SSMs. 
To mention a few of them, we could remember: 

- large enhancement of Ne abundance [138|. 1139] . important because of its contribution uncertainty 
and a weak bond in solar abundances because its abundance is determined rather indirectly [24J; 

- increased element diffusion rates |44l 1140] ; 

- accretion of metal-poor material leading to a 'two-zone' solar model in terms of composition |140|, 

mum]. 

Also solar models including some sort of prescriptions to account for rotation and other dynamical 
effects have been put forward, however their performance is quite poor. 

So far, all attempts of finding a solution to all the manifestations of the solar abundance problem 
have failed. In some cases Y$ can be brougth into agreement with helioseismology, in other cases Rcz 
and the sound speed profile, but a simultaneous solution to all the problems has not yet being found. 

The exceptions are the two obvious one: a) the low-Z solar abundances actually underestimate the 
true metal content of the Sun; b) an increase of radiative opacities by the right amount (15% to 20% 
at the base of the convective zone down to about 3% in the solar core) to compensate for the decrease 
induced by the low-Z abundances. The drawback to this idea is that current state-of-the-art radiative 
opacity calculations differ by only 2 to 3% at the base of the convective envelope, much lower from 
what would be required by low-Z models. 

It has indeed been shown that by increasing the radiative opacity in low-Z SSMs the agreement 
with helioseismology can be restored to match results from high-Z SSMs |142| . Addtionally, the Be 
and 8 B fluxes of a low-Z SSMs with increased opacitiy coincide with those from a high-Z model [143] . 
As good as this may seem, it shows the intrinsic degeneracy between composition and opacities. 

Recently, a novel approach, the Linear Solar Models (LSM), that relate changes in solar observables 
to modifications in the input physics by the calculation of kernels based on SSMs has been developed 
by |144j . LSMgp^j have been applied in particular to the solar abundance problem and the changes 

17 We remark that LSMs offer an efficient way of studying the response of the solar structure to changes in any of the 
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required in the radiative opacity to restore the agreement between low-Z models and helioseismology 
|145j . Quantitively, results similar to those quoted above. 

By using 8 B and now 7 Be as thermometers of the solar core [57], CNO neutrinos represent a unique 
way to break this degeneracy and provide an independent determination of the CNO abundances, 
particularly the C+N abundance in the solar core. Keeping in mind the antagonism between solar 
interior and solar atmosphere models that the solar abundance problem has established, results from 
CNO fluxes will be of the outmost relevance for solar, and by extension stellar, physics. 

14.8.2 The vacuum to matter transition? 

Solar neutrino experiments already measured the two extreme flavor conversion regimes, the vacuum 
term domination and matter term domination. There is no direct experimental evidence of the tran- 
sition from one to the other. In fact, the lower energetic 8 B neutrinos are sensitive to the rise of the 
spectrum from matter combination towards vacuum, but the data (still very uncertain) seems not to 
show it. More data coming from Super-Kamiokande, Borexino and SNO+ experiments will further 
explore the conversion in this regime. 

The precise measurement of low energy neutrinos like pep, exploiting the fact that are more 
energetic than 7 Be neutrinos, will also help to see small solar matter effects in the flavor conversion. 
This matter effects will be more precisely determined by the comparison of pep and pp neutrino 
measurements. In fact, the low energy neutrinos that are better suited to test matter effects are the 
CNO neutrinos. While the CNO neutrinos energy is around the pep neutrinos energy, the former are 
produced at higher temperatures and therefore at higher densities. The larger matter density where 
neutrinos are produced leads to larger matter effects for CNO neutrinos than for pep neutrinos. In 
fact, matter effects produce a significant spectral tilt of CNO neutrinos (~ 10%), which might be a 
good handle to separate the signal for background. 

The determination of the vacuum to matter transition has a significant impact on the determina- 
tion of the solar mass splitting derived by solar data, which adds to the implications of earth matter 
effects measured by comparing the neutrino fluxes during the day and night. The good match of the 
independently determined solar mass splitting by solar neutrino experiments and by reactor exper- 
iments leads to the best test on non-standard neutrino physics to solar neutrinos. There are many 
possibilities but the two scenarios more studied are the addition of new neutral current interactions 
|113| I146| which modify the amplitude of matter effects and therefore shift the effective mass splitting 
and the existence of a sterile neutrino which adds a new state with the appropriate mass splitting 
|114j to produce deviations of the flavor conversion in the 1-3 MeV range. 

14.8.3 What else can we learn from CNO fluxes? 

The most fundamental information the CNO fluxes carry is the most obvious one: that the CNO- 
bicycle operates in stars and it is a viable process for hydrogen fusion. It must not be forgotten 
that neutrinos are the only direct evidence of nuclear reactions being the source of energy in solar 
(stellar) interiors. For the Sun, models predict a marginal contribution to the total energy budget 
from CNO reactions, 0.8% and 0.4% for high-Z and low-Z solar models. However, CNO becomes the 
dominant mode for hydrogen burning in stars with masses right above the solar value. Detection 

physical ingredients entering solar model calculations that does not require the construction of solar models with the 
varied physics. 
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of CNO neutrinos will provide direct evidence that CNO reactions actually take place in nature, as 
originally envisioned by Hans Bethe [147 ] ; it has been a long wait. 

The second important aspect of CNO neutrinos is the information they provide about the abun- 
dance of metals in the solar core. Knowing the abundance of CNO elements in the solar core is 
important by itself. In particular, a 'perfect' measurement of the combined 13 N+ 15 flux translates 
into a determination of the solar C+N abundance with ~ 10% uncertainty [57], and the dominant 
sources of uncertainties are experimental and can be potentially reduced. Assuming we know the 
solar surface abundance of the same elements, i.e. let us forget for the time being about the solar 
abundance problem, we can then put constraints on mixing mechanisms that may have created com- 
position gradients during the evolution of the Sun. SSMs predict that the number density of C+N is 
enhanced in the solar core, at present-day, by ~ 16% with respect to the surface due to the effects 
of microscopic diffusion. And, although helioseismology shows that models with diffusion work much 
better than models without, there is no direct evidence of how efficient diffusion is. In fact, there 
have been suggestions that the standard prescription |148| may be too efficient in the Sun |149j and 
that diffusion rates should be lowered by ~ 15%. Solar CNO neutrinos could provide a test for the 
efficiency of diffusion. 

There are other possibilities that might create a contrast between the solar core and surface com- 
position. Recently, it has been shown that the Sun has a peculiar composition when compared to 'solar 
twins', that is stars almost identical to the Sun in their surface properties }150|. 1151] . The authors 
found that the Sun is enhanced in volatile elements with respect to the solar twins that show no sign 
of harbouring planets by about 20%. In fact, they have associated this fact to the presence of rocky 
planets in the Solar System, where refractory elements are locked, and the occurrence of an accretion 
episode of volatile-enriched material [151] after rocky cores are formed in the protoplanetary disk. If 
this were true, then the Sun would have an envelope that is richer in CNO than its interior. If a 
measurement of the 13 N+ 15 flux would yield as a result a core composition where the abudance of 
C+N would be comparable or less than the surface value, then we would have an extremely exciting 
piece of evidence about the earlier phases of planet formation in the solar system [57] . 
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